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1 Introduction 



This section is devoted to QCD generators, relevant for LEP 2 processes where hadrons may 
be found in the final state: e"'"e~ — > 7*/Z° qq, e"'"e~ W'''W~ qqq'q', e"*'e~ Z°h° — >• 
lyVhh, etc. In fact, almost all interesting processes at LEP 2 give hadronic final states, ensuring 
that QCD generators will remain of vital importance. 

It is instructive to contrast the EW and QCD generator perspectives for LEP 2. In the EW 
physics program, the main emphasis is on four-fermion final states. This is different from LEP 1, 
where the Z° line shape was a major focus of attention |Q. Dedicated four-fermion generators 
are new creations, that have to stand on their own and cannot be tested at LEP 1. Therefore 
there is little sense of continuity with respect to the LEP 1 workshop and subsequent LEP 1 
activities. QCD physics, by contrast, extrapolates logically from LEP 1. New aspects may 
enter, such as colour reconnection, but these are expected to be relatively small perturbations 
on the basic picture (though of importance for precision physics). Therefore the QCD generators 
write-up for the LEP 1 workshop P] is still partly relevant and subsequent LEP 1 experience 
very much so. The high Z° statistics will make LEP 1 a significant testing ground for many 
new QCD physics ideas also in the LEP 2 era. 

It is thus logical to begin this section with an assessment of experience from LEP 1, with 
emphasis on areas where generators are known to have shortcomings. Any improvements for 
LEP 1 will directly benefit LEP 2. This is followed by a comparison of extrapolations to LEP 2 
energies, from which the current range of uncertainty can be estimated. Next comes a survey of 
existing generators, ranging from major programs, with coverage of the full generation chain, 
to shorter pieces of code for specific purposes. Finally, there is a section on standardization 
efforts, to help ease life for users who rely on several generators. 

This report is not a complete description of the topic. However it should provide a convenient 
starting point, with ample references to further relevant literature. 
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2 Experience from LEP 1 



2.1 Event shapes and inclusive distributions 

A large quantitative improvement in the description of event shape and inclusive distributions 
has been made at LEP 1 with respect to the era of PETRA and PEP. This is due mainly 
to the vast amount of high quality data available and the need to achieve good agreement in 
model/data comparisons so as to obtain small systematic errors for the high precision elec- 
troweak measurements. To help facilitate this goal flexible fitting algorithms were developed, 
based on previous work |^. In many cases the dependence of the model's response to its pa- 
rameters is analytically interpolated H, 0, H, |^ . This strategy is flexible, allows easy exchange 
of input distributions but also the simultaneous fitting of very many, 10-15, parameters |^ p!0| . 

Evidently the choice of input distributions used to constrain the model parameters is impor- 
tant. In general a distribution depends on very many parameters, thus the parameters resulting 
from a fit are in general correlated. A survey has been undertaken to determine which distri- 
butions have the highest sensitivity to the individual model parameters |P, 11|. It turns out 



that semi-inclusive spectra are most important, as has been observed before The charged 
particle momentum and transverse momentum spectra strongly constrain the fragmentation 
function or, alternatively, the cluster parameters. However their dependence on the fragmen- 
tation parameters is not exclusive. Inclusive distributions may depend even more strongly on 
Aqcd and/or the parton shower cut-off. In fact, the latter parameter strongly influences the 
high-momentum tail of the momentum spectrum. The 3-jet rate as defined using the Durham 
or JADE algorithms almost only depends on Aqcd- This emphasizes the reliability of as de- 
terminations using this quantity. In contrast, and somewhat surprisingly, the AEEC depends 
strongly on very many model parameters. Measures of the general event topology, e.g. thrust 
and sphericity, depend mainly on Aqcd and only in the 2-jet regime on fragmentation param- 
eters. Shape measures sensitive to radiation out of the event plane, like minor or aplanarity, 
show strong dependence both on fragmentation parameters and on Aqcd- In summary, model 
parameters are best determined by fitting the model to inclusive distributions, jet rates and 
shape distributions simultaneously. 

It appears that the "partonic" phase of the models is best tested by studying the properties 
of jets defined using jet algorithms ||T^. At large resolution parameter i/cut, when dealing with 
few jets or the emission of the "first" hard gluons at large angles, fragmentation effects are 
almost negligible. In contrast at smaller i/cut, where higher jet rates are sizable i.e. when the 
subjet structure described by multiple emission of soft and collinear gluons is important, also 
fragmentation effects are of increasing importance. 



The parton shower models Ariadne |]Tj , HERWIG [^J and Jetset ||T^ describe well the 
general evolution of the individual jet rates with ?/cut, especially the 3-jet rate (see Fig. ^ 
161 and Fig. |] 0). A more detailed 3-jet Dalitz plot study using the ordered normalized jet 



energies Xi {i = 1,2,3) and Z = {x2 — x-i)/V^ unveiled slight differences among the models 
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Figure 1: Differential n-jet rates com- 
pared to Ariadne, HERWIG and Jetset 
PS. 



Figure 2: Differential 2-jet rate compared 
to Ariadne, HERWIG and Jetset PS & 
ME. 



ITf . Ariadne is in perfect agreement with the data, Jetset is slightly below the data in the 
almost 2-jet like case {Z l/"\/3) and slightly above when the lower energetic jets have similar 
energies {Z ^ Q). HERWIG shows a somewhat bigger deviation along the diagonal (i.e. the X3 
direction) of the Dalitz plot. Also the 0{a1) ME option of Jetset is in good agreement with 
the data. It is interesting to observe that the agreement is less good when optimized scales are 
used to achieve a better agreement for the 4-jet rate. 

The discrepancies observed are due either to the different shower evolution strategies used 
or can be traced back to way in which the PS models perform the matching of the first splitting 
to the first order matrix element. In Ariadne this matching is performed naturally, since 
the splitting function is just the lowest order matrix element expression. If no matching is 
performed (a possible option in HERWIG and Jetset) the agreement with the data is poor. 

The 4- and 5-jet rates predicted by Jetset PS (HERWIG) decrease more (less) rapidly 
with i/cut than the data (see Fig. p. At large ?/cut the discrepancy is up to 20% 0. Ariadne 
however is in perfect agreement with the data. 

Clear discrepancies have been observed at PEP and PETRA comparing the Jetset ME 
model to the 4-jet rate. This discrepancy has been resolved by introducing optimized scales 
rS. Today using optimized scales the 4-jet rate is perfectly described by the Jetset ME 
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model 0, ^. However the 5-jet rate predicted by the ME model, as is to be expected, decreases 
far too rapidly and is one order of magnitude below the data at large ?/cut- Recently it has 
been shown that the 5-jet rate is also correctly described [110] when the 0{al) tree-level graphs 



are included in the model ||2^. The scale in this case can be chosen similarly to that for the 
standard O^a"^) case. 

The observations made for the jet rates consistently lead to the following picture if the 
models are compared to event shape distributions: general event shape measures, mainly sensi- 
tive to hard gluon radiation, like thrust, sphericity, Ml^^^^^^^/ ^/s or -Bmax/sum, are reproduced 
extremely well by all PS models |p. The only significant discrepancy is a slight overestimation 
of very spherical events by HERWIG. Observables sensitive to higher order radiation like minor, 
aplanarity, Mj^^/ and -Bmin are consistently overestimated (underestimated) by HERWIG 
(Jetset) for large values of the observables. Due to the normalization of the distributions 
this must also lead to (in general smaller) deviations at intermediate or small values of these 
observables. For example the minor distribution in the case of HERWIG is predicted to be too 
wide. Ariadne is in perfect agreement for most distributions. As Jetset and Ariadne both 
use the Jetset string fragmentation model, it is evident that the discrepancies observed for 
Jetset are due to the parton shower part of the model. 

The general fragmentation part of the models are best tested using inclusive charged particle 
distributions which depend strongly on the interplay between the partonic and fragmentation 
phases of the models. The average charged multiplicity {rich) is the integral of the scaled 
momentum (x) distribution. Both quantities have to be described simultaneously by the models. 
When fitting only to the scaled momentum spectrum, HERWIG predicts {rich) ~ 20.8 close to 
the the very precisely known LEP 1 average {rich) = 20.92 ± 0.24 [^. Ariadne and Jetset 
PS give values that are too small (~ 20.3) and Jetset ME gives too high a multiplicity 
(^ 22.7) il. 



The HERWIG x distribution oscillates slightly around the data distribution. For small x 
it is below, for 0.3 < x < 0.7 it is above (max. 10%) and for larger x again below the 
data. If the multiplicity is constrained to the measured value, the x spectrum is well described 
by the Jetset PS and Ariadne for x < 0.5 but drops 20%-30% below the data for large 
X. This should not to be overinterpreted because experimental smearing is important in this 
momentum range and systematic errors increase. The data so far available from ALEPH and 
DELPHI [^, ^ agree here only within the full experimental error. The Jetset ME result 
also oscillates slightly around the data curve (±5%). 

Thus the multiplicity distribution is described well by Ariadne and Jetset (compare 
Fig. ^). HERWIG predicts a slightly too wide distribution thus overestimating the dispersion of 
the number of hadrons; in HERWIG this is strongly coupled to the number of primary partons. 

The transverse momentum in the event plane, p±m, is strongly sensitive to hard gluon 
radiation and almost correctly described by all models. Only the large p^in tail is slightly 
underestimated. The predicted p±out distribution for p±out > 0.8 GeV falls off more rapidly 
than the data in all models and at large p±out is ~ 30% below the data! This fact is shown in 
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Figure 3: Distribution of p±out with re- 
spect to the sphericity axis compared to 
Ariadne, HERWIG and Jetset PS & 

ME \m. 



Figure 4: Multiphcity distribution 
compared to the DELPHI tuning of ARI- 
ADNE, HERWIG and Jetset PS & ME 



figure y which also compares the data of ALEPH and DELPHI [^] to depict the precision 
of the experimental data. The large p±ont tail is mainly due to gluon radiation. This failure 
of the shower models is presumably due to missing large angle contributions in the basic LLA 
used by the models. A matching of the second order matrix element and the LLA shower 
formalism should lead to an improved description similar to that of the matched NLLA and 
calculations used in as determinations |^ 



For the ME model the situation can be 



improved by including higher order terms as has been shown recently by OPAL [19 



2.2 Particle composition and spectra 

Experimental studies of the spectra and composition of particles in hadronic jets provide an 
unique way to understand the fragmentation of quarks and gluons into hadrons. Thanks to the 
excellent performance of the detectors and high statistics available, very careful work by all four 
LEP experiments has given us a very complete picture of the production of identified particles 
from e^e~ annihilation. All states in the SU (3) pseudoscalar and vector meson nonets, except 
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the p"*", and at least one state per isospin multiplet in the baryon octet and decuplet, plus 
the scalar fo(980) and the tensors f2(1270), K2(1430) and f2(1525) |^ have been measured. 
The average production rates per hadronic Z event, together with the predictions from the 
tuned Jetset 7.4 and HERWIG 5.8, are listed in table |l| The measurements are 

in good agreement between experiments for all mesons and octet baryons. However for the 
decuplet baryons there are still discrepancies between experiments, reflecting difficulties in 
the measurements. In particular, the A"*""^ signal is difficult to measure because of its large 
width and the large combinatorial background in combinations of 7r"'"p. The Q~ rate seems 
to be established around the value expected from Jetset, contrary to the old claims of an 
anomalously high production rate. 

Particle rates could depend on many things, such as flavor content, spin, mass, phase space, 
hadron wave functions, Bose-Einstein interference and other collective effects. The two most 
frequently used models HERWIG and Jetset use different ways to account for the particle 
production rates. In the Lund/JETSET approach (similarly to the old Field & Feynman model 
the production rate of a specific hadron type depends principally on its flavor content 



and spin. One can also use essentially pure phase space as in the case of the HERWIG cluster 
fragmentation approach. 

Studies of general features of particle production, such as the strangeness suppression factor 
s/u or the fraction of mesons produced in spin-1 states, V/ (V + P), or in orbitally excited states 
provide useful information about the main production mechanisms. The (one dimensional) 
string model suggests the production of orbitally excited states is small ~10% |^ whilst 
V/{V + P) = 3/4 is expected from simple spin counting. 



Measurements of the fo(980), f2(1270), K;(1430)° [0 and f^(1525) [|| as well as of D 



and B** mesons indicate that orbitally excited states, most of which so far were not included in 
HERWIG, Jetset and other models, are copiously produced (~ 30% of the primary hadrons). 
Thus a quite large fraction of the observed stable particles come from decays of these numerous 
states. As a result, the V/{V + P) ratio can differ signiflcantly from that when no orbitally 
excited states are considered. From a global tuning, where the orbitally excited meson states 

are included, a value of V/{V + P) ~ 0.4 0.6 is obtained for light mesons |2D[. This low 

value ofV/(y + P) could be explained by mass differences between the vector and pseudoscalar 
mesons, i.e. by the relatively larger binding energy of pseudoscalar mesons [P^. The measured 
ratio of V/{V + P) = 0.75 ± 0.04 [^] for B mesons agrees well with the expected value of 3/4. 



However for D mesons the much lower value 0.46 ± 0.06 Pljl is still not understood. 



In the string fragmentation model, one expects the strangeness suppression factor s/u to 
be around 0.3 using the typical values of (constituent) quark masses. This parameter can be 
measured from the production rates of strange compared with non-strange mesons and from 
the momentum spectrum of strange mesons. The results, which are summarized in table |^, 
are very consistent with the expectation^. It is interesting to see that s/u determined from 



^However, neutrino experiments at lower energies |32 and recently both ZEUS |Q and HI |Q require a 



lower value of about 0.2 for s/u. More careful studies in this area are needed in the future. 
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Table 1: Average particle production rates in hadronic Z decays (excluding charge conjugates 
and antiparticles) , compared to the predictions of Jetset and HERWIG. A * indicates that 
the predicted rate differs from measurement by more than three standard deviations. 
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heavy mesons agrees well with the values obtained from light mesons. This suggests that the 
strangeness suppression occurs at the quark level. 



Technique 



Results 



References 



K*° K'^ 20 2l 

at high momentum 
at high momentum 



0.29 ±0.03 

0.35 ± 0.07(stat) 
0.25 ± 0.03(stat) 



table |1| |5[ 



|35|, HI 



|35|, |7| 



K*^ momentum spectrum 


0.285 ±0.035 


|38| 


K*^ momentum spectrum 


0.30 ± 0.02(stat) 


|39| 
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0.32 ±0.08 

0.28 ±0.11 
~ 0.3 
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Table 2: Measurements of s/u at LEP 



The relative rates of the decuplet baryons as well as those of the and provide a 
clean test of models, since they are less affected by resonance decays. From table |1], we obtain 
the following ratios: 



Ratio 


Measured 


Jetset 


HERWIG 




0.29 ±0.04 


0.44 


0.84 




0.36 ±0.07 


0.42 


0.93 


E(1385)+/A++ 


0.17±0.09 


0.11 


0.21 


S(1530)7S(1385)+ 


0.28 ±0.07 


0.38 


0.77 


fi-/H(1530)° 


0.26 ±0.09 


0.28 


0.31 



One obtains from the above ratios a suppression factor of about 0.28 per s quark for baryons 
(0.24 if only the decuplet baryons are considered). This is similar to the value obtained for 
mesons, suggesting that the additional suppression for diquarks might be small. 

After being tuned to LEP 1 data, HERWIG and Jetset^ describe well the measured rates 
in the meson sector. There is a fairly good agreement in the baryon octet, except that the 
proton rate is slightly underestimated and the S~ rate is overestimated by about a factor of 

^New parameters have to be introduced, as attempted in to treat the quark type dependent production 
probabihties for pseudoscalar, vector and orbitally excited mesons. 
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two by HERWIG. In the baryon decuplet Jetset predictions are consistent with the data while 
the predictions of HERWIG differ from the data in most of the cases. Differences in the ratios 
of the baryon rates between HERWIG and data, as shown above, can not be solved simply 
by tuning the cluster fragmentation parameters, indicating the need for real dynamics beyond 
phase space and spin counting. 

Although in general Jetset describes the measured rates better than HERWIG, it contains 
a large number of free parameters. As a result, it has little predictive power. The UCLA model 



44| , a variant of Jetset with less parameters, does a good job in many cases but has problems 
in the baryon decuplet. Also the copiously produced orbitally excited mesons so far are not 
included in the UCLA model. In |^ an interesting regularity in production rates is shown 
for all particles (except pions) belonging to the pseudoscalar and vector meson nonets and the 
baryon octet and decuplet. The particle multiplicity can be described by a simple exponential 
fall off in mass squared and 2 J + 1 spin counting factors. This regularity seems to be energy 



independent and has recently been established similarly also in pp interactions [^. However 
it is necessary to use generalized isospin multiplets and to not seperate the contributions from 
resonance decays. Recently a new approach |^ has been proposed which uses only three free 
parameters but reproduces the measured rates quite well. The basic assumption used is that 
hadrons reach complete thermal and chemical equilibrium, in contrast to the general belief that 
e"'"e~ — >-hadrons is a rapidly expanding process and during fragmentation the hadronic density 
is rather low. More tests are needed to check this thermodynamic approach. 

Since all fragmentation models contain a number of parameters which can be tuned accord- 
ing to data (more dramatic in the case of Jetset), measurements of production rates do not 
provide a high discriminating power among different models. A more effective method is to 
look at baryon correlations. In Jetset the major source of baryon production is the creation 
of a diquark-antidiquark pair within the fragmentation. The baryon-antibaryon (BB) rate is 
much higher than the BB and BB rate (see table |I]) and BB are more likely to occur close in 
phase space than BB or BB. Correlations between BB can be reduced by the popcorn mech- 
anism, allowing a meson to be created in between a BB pair. As can be seen from table 
HERWIG overestimates the AA rate (note that the prediction for the A is quite good), while 



Jetset with popcorn describes the data well. It has been shown in that BB correlations, 
for example in rapidity, are overestimated by HERWIG, whilst Jetset with a high probabil- 
ity of the popcorn occurrence (~ 80%) reproduces the data well. A more impressive test is 
to study the angle between the baryon and the event axis in the BB rest frame. The string 
model predicts that baryon production is preferentially lined up along the event axis, while the 
cluster model predicts an isotropic distribution. Data |^ clearly favor the string model. Also 
measurements of baryon and antibaryon production in quark and antiquark jets with polarized 
beams by SLD |^ and jet charge studies support the string model but disfavor the idea 
of isotropic cluster decays. 

Identified particle spectra have been studied as function of both the scaled momentum 
= p/Pbeam and the variable = log ^. In general all models describe the data fairly well, 
with few discrepancies remaining: 
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Data show a harder momentum spectrum for the t] produced in gluon jets 
K='= momentum spectra predicted by the models are too soft |3^, 0]. This might be 
caused by wrong branching fractions of b hadrons in the models. 
Momentum spectra of light quark baryons predicted by the models are too hard 

1371 m 



3e, 



m 52 



This indicates a different production mechanism for baryons than for 
mesons. Partly it may also be due to missing orbitally excited baryon states in the 
models. However, so far these states have not been observed in e"'"e~ annihilation. 
The heavy quark fragmentation function has been measured at LEP 1 mainly using D^*^ 
reconstruction in the c-quark case and using high-p^ lepton spectra D*— lepton com- 
binations and exclusive and inclusive b- reconstruction |5^ in the case of b-quark 
fragmentation. The D-meson distributions are obscured by contributions from b-hadron de- 
cays. Today the (experimentally involved) inclusive b-hadron reconstruction yields the best 
statistical precision. It allows for the first time (besides a precise determination of the average 
b-hadron energy (xe)) a decisive comparison to different fragmentation models. This, so far 
incomplete comparison, gives best agreement for LLA based parton shower models (Ariadne, 
Jetset and HERWIG) combined with Peterson fragmentation [^]. The HERWIG cluster frag- 
mentation as well as the Lund-symmetric and the modified Lund-Bowler ansatz give less sat- 
isfactory results. In the case of Jetset ME with Peterson fragmentation a too narrow energy 
distribution indicates the lack of soft gluon emission. 



2.3 Differences between q and g jets 

In QCD, the gluon is associated with a color charge Ca = 3 and the quark with a charge 
Cf = 4/3. The larger color charge of the gluon means that it is more likely to radiate an 
additional gluon than a quark, leading to differences in the expected properties of quark- and 
gluon-induced jets. For quark and gluon jets produced with the same energy and under the 
same conditions, gluon jets are expected to have a larger mean particle multiplicity than quark 
jets [^. The larger multiplicity of the gluon jet implies that its particle energy spectrum. 



known as the fragmentation function, is softer. A related prediction is that the mean opening 
angle of particles in a gluon jet is larger than in a quark jet [Q: thus the gluon jets are 
broader. Much experimental effort has been invested in an attempt to observe these predicted 
differences (for a recent compilation, see and references therein). Before LEP 1, there 



were experimental indications that gluon jets were indeed broader than quark jets, based on 
measurements of the mean transverse momentum of particles in a jet with respect to the jet 
axis, or similar variables. However, contradictory results were published concerning differences 
between the quark and gluon jet fragmentation functions, while no evidence was found for 
a multiplicity difference between the two jet types. In general, it proved difficult to obtain 
conclusive results on quark-gluon jet differences at facilities before LEP 1 either because biases 
were introduced by assuming the gluon jets to be the lowest energy jets in e+e~ three-jet events 
or else because there was no event-by-event identification of gluon jets with a resulting lack of 
sensitivity. 
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Due to large event statistics and good detector capabilities, the LEP experiments have been 
able to settle the experimental question of quark and gluon jet differences 33]. Three 



aspects of the LEP 1 studies allow this success. (1) Symmetric events were selected in which 
the quark and gluon jets being compared had the same energy and angles relative to the other 
jets, allowing a direct, model independent comparison of the jet properties. (2) The quark jets 
were tagged, leading to identification of the gluon jets with better than 90% purity through 
anti-tagging. (3) The anti-tagged gluon jet data were combined algebraically with the quark 
and gluon jet data from the untagged, symmetric events, leading to separated quark and gluon 
jet measurements with essentially no biases except from the jet definition. In the first LEP 1 
studies, the quark jet samples were the natural ones for T? decay, given by the Z° coupling 
strength to the individual flavors, corresponding to roughly 20% d, u, s, c and b quarks. In 
later studies, b quark jets and uds quark jets were explicitly selected to compare to gluon 
jets [§4 . 



These studies resulted in a confirmation of the qualitative differences between quark and 
gluon jets given above. Selecting 24 GeV jets in a so-called "Y" symmetric event topology, it 
was shown that gluon jets were 60-80% broader than quark jets as measured by the full width 
at half maximum of the differential energy and multiplicity profiles The fragmentation 

function of the gluon jet was observed to be much softer than that of the quark jet. The mean 
charged particle multiplicity of gluon jets was found to exceed that of quark jets by 20-25%. 
Besides the Y events, DELPHI studied 30 GeV jets from three-fold symmetric "Mercedes" 



events and obtained similar results. The comparison of the fragmentation function of quark 
and gluon jets in Y and Mercedes events shows the expected stronger energy dependence for 
gluon jets. Extensive comparisons of Monte Carlo predictions to the quark and gluon jet data 
are presented in and Ariadne, HERWIG and Jetset were found to be in good 



agreement with the measurements. The COJETS agreement was somewhat less good. 



ALEPH 1^ extended these studies by including a measurement of sub-jet multiplicities p7| . 
For small values of the sub-jet resolution scale, (defined using the jet finder), the ratio 
of the gluon to quark jet mean sub-jet multiplicity was found to be similar to the hadron level 
value of about 1.2 discussed above. After subtracting one from the mean sub-jet multiplicities 
to account for the contributions of the initiating quarks and gluons, the sub-jet multiplicity 
ratio of gluon to quark jets was observed to reach a much larger value of about 2.0 as 
approached the resolution scale y\ at which the jets were defined. The explanation for this is 
that the mean sub-jet multiplicity of the quark jets approaches unity slightly before that of the 
gluon jets as y^ y\. Ariadne, HERWIG, Jetset and NLLjet were all found to reproduce 
the measurement. 

Beyond these studies based on symmetric events, ALEPH and DELPHI have examined 
quark and gluon jet properties in non-symmetric three-jet event configurations. The DEL- 
PHI approach [^| is to identify gluon jets in three-jet events using anti-tagging methods as 



mentioned above. The gluon jet properties were compared to those of quark jets with similar 
energies found in radiative QED qq7 events. The qualitative differences discussed above be- 
tween quark and gluon jets were observed to be present for jet energies between 5 and 40 GeV 
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and were well reproduced by Jetset. ALEPH [^] introduced a new method to study the 
multiplicity difference between quarlc and gluon jets in three-jet events, by examining the mean 
charged particle multiplicity of the entire event as a function of the energies and opening angles 
of the jets in the event. Assuming each event to be composed of a gluon jet and two quark jets, 
and that every particle in an event could be associated with one of these jets, a fit was made to 
extract a value for the ratio of the mean charged particle multiplicity values of gluon to quark 
jets, Teh- The result for all jet energies and event topologies was Teh = 1.48. The fit results 
were found to agree well with those from the symmetric Y analyses when they were restricted 
to that geometric situation. 

Thus the basic differences expected between quark and gluon jets — a larger mean multi- 
plicity, a softer fragmentation function and a larger angular width of gluon relative to quark jets 
— are now all well established by the LEP 1 experiments. The QCD models are in good overall 
agreement with the measured differences. Future effort in this field at LEP 1 will probably in- 
clude studies of differences in the identified particle rates in gluon and quark jets, differences in 
particle correlation phenomena and attempts to reduce the reliance of the analysis method on 
the jet definition (as the ALEPH study discussed above attempts to do). Already, L3 has 
presented results which indicate an enhanced r] meson production rate in gluon jets compared 



to the rates predicted by HERWIG and Jetset |50]. This suggests that the models for gluon 



jets may need to be modified to allow for an enhanced production of isosinglet mesons [3£]. 



2.4 Coherence 

Gluon radiation in the parton shower should be coherent. However, gluon interference only 
becomes apparent when one goes beyond the Leading Log Approximation (LLA). A number of 
such effects are found in the next simplest approximation, the Modified LLA (MLLA) [^]. Due 
to the non-abelian nature of QCD, the overall result of this interference is "angular ordering" 
of the gluon radiation [|7l|], which constrains the angles between the radiator and the radiated 
gluon to decrease as the evolution proceeds to lower scales. 

In parton-shower Monte Carlos gluon interference is either: imposed as an a posteriori 
constraint on gluon opening angles as in Jetset [T^ ; built into the choice of evolution variable 



as in HERWIG [T^; or neglected in independent fragmentation models such as Co jets |72 
Ariadne |TB[, on the other hand, employs a formulation based on a cascade of qq, qg and gg 
dipoles which naturally incorporates interference phenomena. In Jetset the angular-ordering 
constraint can be turned off. By comparing Jetset with and without angular ordering one 
can obtain an idea of the importance of the effect. 

Some consequences of gluon interference have been calculated directly in perturbative QCD 
as well as by Monte Carlo. Such calculations apply, strictly speaking, only to partons. Com- 
parison with data relies on the additional assumption of Local Parton Hadron Duality (LPHD) 



[[73i [74| , which posits that many distributions of hadrons rather closely follow the correspond- 



ing parton distribution, with non-perturbative effects affecting mainly the normalization rather 
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than the shape of the distributions. However, we shall not emphasize such calculations here, 
since our main purpose is to evaluate the adequacy of current Monte Carlo programs. 



The first effect to be explained |7§] as a consequence of gluon interference was the so-called 
string effect; first predicted using (non-perturbative) string fragmentation phenomenology 
and later discovered by the JADE experiment . In terms of gluon interference it is explained 
as a purely perturbative effect at the parton level. The string effect has been extensively studied, 
most recently by DELPHI L3 and OPAL ||79[. These analyses have compared qqg 
and qq7 events taking care to have samples of comparable kinematic configurations. The string 
effect appears as a smaller particle flow in the region between the quark jets in qqg than in qq7 



events. ALEPH ||8^ instead compared the particle flow between the quarks with that between 
quark and gluon. The string effect is found to be rather well reproduced by the coherent Monte 
Carlo models but not by the incoherent ones. However this success is not entirely due to the 
coherence at parton level; the non-perturbative modelling in the programs also contributes. 

It is also worth mentioning that evidence of gluon interference is also seen in pp interactions 
at the Tevatron. Using events with 3 high-p^ jets CDF examined the differences in rapidity 
and in azimuthal angle between quark and gluon jets ||8l|. HERWIG, which incorporates co- 



herence in both space-like and time-like showers, reproduced the data well. Pythia/Jetset, 
with coherence only in time-like showers did less well, although it improved when modified to 
partially incorporate coherence in space-like showers. Isajet, with no coherence, performed 
poorly. 



As is well known |7y], gluon interference leads to suppression of soft gluons in the shower, 
which in turn should lead to a suppression of soft hadrons. The distribution of = — In Xp = 
— Inp/pjet is expected to have a roughly Gaussian shape and its peak position, should 
increase with ^/s. The dependence of ^* on ^/s is strikingly different in the MLLA from that in 
the LLA. Assuming LPHD, C,* is expected to show similar behaviour. Many comparisons have 
been made, for many types of particle, using data from PETRA/PEP, TRISTAN, and LEP; 
they support the form predicted by MLLA and clearly reject the LLA form. 



From MLLA-I-LPHD it is expected |Q that ^* decrease with the mass of the hadron. This is 
indeed found to be the case with ^* being approximately proportional to — InMhadmn- However, 
the proportionality constant is quite different for mesons and baryons This difference is 



due to decays. When the distributions are corrected for decays using Jetset, 

the ^* values of mesons and baryons are found to lie on a universal curve |^]. The conclusion 
is clearly that we must be cautious about the interpretation of LPHD, in particular with the 
inclusion of decays. 

The y/s dependence of ^* is support for MLLA, but says little about the quality of the 
Monte Carlo programs, since they are retuned at each value of ^/s. However, accepting the 
validity of MLLA, the improvement seen in the previous paragraph supports the description of 
non-perturbative hadronization in the model. 

The angular ordering resulting from gluon interference effectively moves the radiated gluons 
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closer to the jet axis. The size of the effect depends on the colour charge of the radiator and 
on the initial configuration of the event (i.e. 2 or 3 jets, there being interference effects in 
the interjet region for 3-jet events as seen in the string effect). The total number of (sub-)jets 
found in an event has been calculated in the Next-LLA (NLLA) as a function of the jet 



resolution parameter, y^ut = Ho-, for 2- and 3-jet events classified using Hcnt = Vi > Vo- Of 
course, if yi becomes too small non-perturbative processes become important and the calcula- 
tion breaks down. The perturbative and non-perturbative regions are rather clearly separated 
and the sub-jet multiplicities thus provide a test not only of perturbative QCD calculations 
and their incorporation into Monte Carlo programs, but also of the non-perturbative models in 
the programs. 

Sub-jet multiplicities have been studied [0, ^ at LEP 1. Quite good qualitative agreement 
is found between the data and the NLLA calculations in the perturbative region while a simple 
0{as) calculation clearly disagrees. Of the Monte Carlo programs, Ariadne does quite well; 
HERWIG 5.5 and Jetset 6.3 perform somewhat less well; and the incoherent model Co JETS 
gives the worst agreement. Both versions 6.12 and 6.23 of Cojets disagree in the perturbative 
region while only 6.23 disagrees in the non-perturbative region. Jetset was compared ||85| using 



various combinations of fragmentation and parton shower schemes. Incoherent parton showers 
resulted in poor agreement in both perturbative and non-perturbative regions independently of 
the fragmentation scheme. Coherent showers gave much better agreement in the perturbative 
region. In the non-perturbative region agreement was poor for independent fragmentation 
whilst good for string fragmentation. 



The MLLA predicts |86| a suppression of gluon emission within a cone of angle 6 < 6q = 
Mq/Eq about the quark direction in a parton shower. This should lead to a lower primary 
multiplicity for heavy quark events. However, the total multiplicity is higher because of the 
high multiplicity of heavy flavour decays. The difference in multiplicity between heavy and 
light quark events is predicted to be independent of y/s, contrary to the naive expectation that 
the difference would decrease as the quark mass difference becomes smaller compared to the 
total energy. Results from PEP/PETRA, TRISTAN, and LEP/SLC agree reasonably well with 
the MLLA value, both for charm and beauty, particularly when the recent work of Petrov and 
Kisselev is taken into account. The results |]88|, |89|, 190] at ^/s = Mz agree reasonably well 



with the predictions of Jetset, with the possible exception of light (uds) quarks. 

Given the appearance of angular ordering in MLLA, the effects of gluon interference should 
be apparent in angular correlations. Assuming LPHD, the correlations should persist in the 
hadrons. Besides the angular ordering in the polar angle, also the azimuthal angular distri- 
bution is affected by gluon interference. OPAL ||9l]], has studied two-particle correlations in 



the azimuthal angle within restricted rapidity intervals. To avoid defining a jet axis they, and 
more recently ALEPH [^, have also studied such correlations using the Energy-Multiplicity- 
Multiplicity Correlation (EMMC) Taking in turn each track's direction as an axis the cor- 
relation calculated, the EMMC is the average of these correlations weighted by the axis defining 
track's energy. The EMMC has been calculated analytically in leading and next-to-leading 



94 1 order; the corrections are large. LPHD must be assumed to apply these calculations to 
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those calculated from hadrons. Nevertheless, qualitative agreement is obtained for > 7r/2, 
where Monte Carlo models show hadronization to be relatively unimportant. Agreement with 
the data is even better for Monte Carlo models which incorporate gluon interference. Models 
not incorporating this interference fail to describe the data. 

ALEPH [p2| , |95[| and L3 |9^, ^ have studied two-particle angular correlations in the full 



spatial angle using the Asymmetry in the Particle-Particle Correlation (APPC). In addition, 
L3 has studied the Asymetry in the Energy-Energy Correlation AEEC. The APPC is defined 
in analogy to the well-known AEEC by simply removing the energy weighting. This results in 
a correlation which is sensitive to all branchings of the shower, whereas the AEEC is primarily 
sensitive to the earliest branchings. The APPC is less sensitive to systematics in the correction 
for detector effects. On the other hand, the energy weighting makes the AEEC less sensitive 
to the Bose- Einstein effect. The use of the asymmetry serves to cancel some of the correlations 
arising from other effects as well as some detector effects and Monte Carlo uncertainties. 

These correlations have been compared with Monte Carlo models. The conclusion is that the 
models containing gluon interference agree much better with the data than do the incoherent 
models. However neither version of NLLjet can be said to agree well. 

All of these studies favour the Monte Carlo models Ariadne, HERWIG and Jetset, which 
incorporate the gluon interference expected in MLLA. In general the agreement of data with 
these models is quite good. On the other hand, models that do not incorporate gluon interfer- 
ence, such as CoJETS and incoherent Jetset do not in general agree well with the data. Both 
coherent and incoherent versions of NLLjET have been found also not to agree well with data. 



2.5 Prompt photons 

The principal source of observable prompt photons in hadronic decays of the Z (i.e. those with 
energies greater than a few GeV) is final state radiation (FSR) emitted at an early stage in 
the parton evolution process initiated by the primary quark-ant iquark pair. To reduce large 



backgrounds from non-prompt sources, the first measurements reported by OPAL ||98[ and 
followed later by the other LEP experiments selected events with photons well isolated from the 
hadrons by a geometrical cone followed by a 2-step jet reconstruction process. In this procedure, 
the candidate photon is first removed from the event and all the hadrons reconstructed into 
jets. Then, the photon is replaced and its isolation from the jets tested in a second application 
of the clustering algorithm. It was soon realized that the cross sections are substantially less 
than those predicted from fractionally charged fermion pairs due to the influence of gluons. 
Thus, the measurement of prompt photons has become a sensitive test of the predictions of 
both perturbative QCD matrix element calculations and the Monte Carlo shower models free 
from the direct effects of fragmentation. 

After tuning the parameters of these models in recent versions, namely Ariadne 4.2, 
HERWIG 5.4 and Jetset 7.3 to the properties of the hadrons observed in non-FSR events. 
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there is no freedom to adjust the photon emission parameters with the exception of the infra-red 
cut-off. In the following reported analyses, these cut-off values are chosen to be similar to those 
employed to terminate the parton evolution, but in any case do not significantly influence the 
isolated hard photon rates. 



All the published high statistics analyses from ALEPH, L3 and OPAL [^] show that both 



Ariadne and HERWIG give acceptable descriptions of the total and individual n-jet -|- 7 
cross sections as a function of the jet resolution parameter, ?/cut (JADE EO), as well as the 
distributions in p± and fractional energy of the photon. A more critical test to differentiate 
between these two models is based on their predictions for the rate of low energy photons (< 15 
GeV) at large angles (> 75°) to the event thrust axis, where the evolution scale ordering used 
in HERWIG predicts a larger cross section than Ariadne | |10CI| | . Preliminary data from ALEPH 



indicate that Ariadne gives the better description but more statistics are needed. However, 
the above published results show that Jetset is less satisfactory predicting cross sections that 
are 20-30% low (Scr's). ALEPH showed that this can be improved by either switching off the 
0{as) matching or by keeping as constant indicating that virtuality as the scale controlling 



the parton evolution is not the best choice. More recently, DELPHI has also shown ||101|| that 
their data are in excess of Jetset by 18±7% in the low energy region of the photon spectrum 
below 15 GeV. After clustering the hadrons with the Durham {k±) algorithm in a similar 2-step 
procedure, their respective jet rates above Hcut ^ 0.01 are in reasonable agreement with Jetset. 
The excess is largely eliminated since most low energy photons are no longer isolated when the 
clustering algorithm is applied a second time. This appears to be a different conclusion from the 
other experiments. However, careful examination shows that the discrepancy between Jetset 
and data for ALEPH and L3 are largely at low i/cut in the total cross section where the use of 
different algorithms for jet-finding makes comparison difficult with DELPHI. It should be noted 
also that DELPHI compare with Jetset at hadron level before fragmentation corrections are 
applied. 

The 2-step analysis procedure to select isolated photon events does not prevent a significant 
number of non-isolated hard photons from contaminating the 7 -|- 1-jet event topology. Each 
of these photons remain within the hadron jet formed from the remnant of the radiating quark 
and are better separated from the isolated radiation by a "democratic" analysis | |102| | . Here, 



the prompt photon candidate is not removed from the event and thus becomes a member of a 
hadron jet with fractional energy of its total energy. The true isolated component is now 
concentrated at = 1 broadened downwards in by hadronization effects to overlap with the 
high energy tail of the coUinear quark fragmentation component. For the 7 -|- 1-jet (ie: 2-jet) 
cross section, this is well separated from the fragmentation tail when z^ > 0.95. Fig. ^ shows the 
comparison as a function of ?/cut (Durham EO scheme) between the data measured by ALEPH 
and the predictions of Ariadne and Jetset for this isolated component. The continuous 
curve is a prediction of a leading order calculation dominated by perturbative terms which are 
derived from a pure QED calculation. HERWIG (not shown) is in close agreement with the 
data. Jetset falls well below the data in this case showing that its treatment of radiation 
as independent emission from either quark at the first branching is quite inadequate. This 
discrepancy diminishes as the jet multiplicity increases. There is more satisfactory agreement 



18 



Figure 5: Integrated 2-jet rate above = 
0.95 function of ycut, compared with 
Ariadne, Jetset and a QCD calcula- 
tion. 



in the fragmentation region below Zj = 0.95. 

Overall, the conclusion is that HERWIG gives the best description of all prompt photon 
data at the Z closely followed by Ariadne. 

In this review it is appropriate to mention the difficulties faced in determining the non- 
prompt photon background coming from hadrons decaying into 7's (mainly 7r°). The isolation 
and energy cuts applied to the prompt photon candidates in the 2-step analyses are insufficient 
to eliminate this background entirely even when the full granularity of the electromagnetic 
calorimeters is exploited to recognize single from multiple 7 showers. Hence, an irreducible 
non-prompt component must be subtracted statistically using QCD models or inferred from 
other data. However, the selection cuts applied choose a region of phase space that is not 
well understood in these models as they correspond to tails in the fragmentation process which 
cannot be tuned precisely. 

The early analyses made at LEP 1 showed a clear discrepancy in the hadronic background 
yield predicted by the HERWIG and Jetset models [[L04|] . The magnitude of the differences 
depends strongly on the isolation and energy cuts. A substantial effort has been made to 
quantify these discrepancies in detail, most recently by L3 [|105|| . They are able to reconstruct 
well resolved vr^s and r^s from two identified photons isolated by a geometric cone in which no 
other particles are found with energies above 50 MeV. Jetset reproduces the observed rate 
of 7r°s and r]s with energies above 3 GeV for 10° isolation, but significantly underestimates the 
rate for 25° isolation. This study was restricted to 8 GeV maximum energy where the direct 
meson reconstruction procedure is efficient, but has been extended to 45 GeV using a neural 
network. The observed background rate of non-prompt photons is about a factor 2 larger than 
the predicted rate over the full energy range and the discrepancy increases with tighter isolation 
cuts. HERWIG tends to give a slightly better prediction but still underestimates the rate. 



In other studies at the Z of the non-prompt photon background both ALEPH |9^] and 
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DELPHI 101 



have reported that Jetset underestimates the isolated 7r° yields but only in 
the lower part of the energy spectrum below 20 GeV. In these analyses, the limit allowed for 
the maximum particle energy accompanying the photon in the cone is set to 500 MeV. They 
are not inconsistent with the L3 results but instead demonstrate that the comparison with the 
generators is sensitively dependent on the isolation parameters. In the alternative "democratic" 
analysis without isolation cones of ALEPH ||102|| some activity is allowed in the vicinity of the 
7 which results in a better description by Jetset of the region of phase-space considered for 
the fragmentation. 



2.6 Bose— Einstein effects 

Most of the Bose-Einstein interference studies at LEP 1 have concentrated on two-particle 
correlations between identical charged pions ||106|| using the quantity 



m) = ''^^L (1) 

Here, P2{M) is the two-particle correlation function, usually given as a function of Q, = 
— 4m^, and pi ® Pi{M) is a reference sample. This sample should resemble p2(M) except 
for the Bose-Einstein correlations being studied. 

Two choices for the reference sample are made, unlike-sign pion pairs or uncorrelated pairs 
from track mixing. Both alternatives have disadvantages. Unlike-sign pion pairs suffer from 
correlations due to resonances not present in like-sign pion pairs and the contribution of reso- 
nances with poorly known rates, especially rj and rj' at low Q. Furthermore residual effects of 
Bose-Einstein interference may also be visible in the unlike-sign pairs (see below). The track 
mixing has the disadvantage that correlations, other than from Bose-Einstein interference, are 
missing. In addition cuts to suppress gluon radiation must be applied. For both methods the 
systematic uncertainties are reduced using the double ratio R'^^^'^{M)/R^'^{M). Additional 
corrections for background, e.g.. Coulomb interactions are applied. 

Assuming a spherical and Gaussian source the enhancement at low Q is parameterized as 
R{M) ~ 1 -|- Aexp(— r^Q^). The chaoticity parameter A is expected to vary between and 1, 
and is extracted from data in the range from 0.4 to 1.5; the radius r of the source is measured 
to be 0.4 fm to Ifm. In Fig. ^ the background-corrected measurements are displayed for the 
mesons tt^, 7r°, K^, and K°. 

Only identical mesons, that are prompt, i.e. do not originate from long-lived resonances, 
can contribute to the enhancement at low Q. It has been pointed out that the measured value 
of A is about the maximum you could expect from direct pairs or even higher [ |107| ]. 

In more recent analyses the fraction f{Q) of direct pions as a function of Q has been 
parameterized using Monte Carlo and included in the fit. For example DELPHI uses f{Q) = 
0.17 + 0.26Q - 0.12Q2, obtained from Jetset, to fit A and r for charged pions: R{M) ~ 
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Figure 6: Chaoticity parameter A versus radius r measured at LEP 1. Measured values are 
corrected for background with statistical (solid line) and total errors (dots) shown. The arrows 
indicate the changes, when corrected for non-prompt meson-pairs estimated with HERWIG or 
Jetset, when it is calculated by the experiment ||106|, |108|, |109||. 



1 + Xf{Q) exp(— r^Q^) [|106|| . Whilst the change in the radius is small, A is changed by a factor 
3. A bigger change is reported by L3 for 7r° — tt" correlations ||108| | . The corrections are very 



sensitive to the model used. The corrections for non-prompt mesons are indicated by arrows 
in Fig. ^. The kaons have higher chaoticity values than pions before correction |l(J9|. Only 



DELPHI has estimated the corrections for non-prompt kaons. The correction for kaons from 
c- and b-decay increases A by ~ 25 to 30%. 

Three-particle correlations have been studied by DELPHI. Whilst Jetset without Bose- 
Einstein correlations fails to describe the data, Jetset with Bose-Einstein correlations enabled 
gives a fair description of unlike-sign triplets; the shape is reproduced, but the magnitude is 
too small ||110| . 



Bose-Einstein correlation affect the unlike-sign spectra as well. In the invariant mass distri- 
bution of pions the p° meson appears shifted towards lower masses [111]. In the framework of 



the model this can be interpreted as coming from Bose-Einstein correlations between like-sign 
pion pairs, which induces correlations between unlike-sign combinations, for example seen as a 
distortion of the p° line shape. OPAL finds nice agreement between data and Jetset including 
Bose-Einstein correlations, when the chaoticity parameter is set to 2.5. This value of A was 
obtained with a fit to the ratio R{M). ALEPH agrees with this observation and extracts a p° 
rate with A and r as free parameters. The value of A = 2.1 is compatible with OPAL in view 
of the different r/' rate and choice of the coherence time parameter x- {x gives the minimum 
width of resonances whose daughters contribute to the Bose-Einstein enhancement). DELPHI, 
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which observes a shift of the uses its A value extracted from the Bose-Einstein analysis, 
after correction, for the p° analysis. Also with this parameter choice the agreement of data and 
model mass spectra is satisfactory ||112| . 



Concerns have to be raised about the implementation of Bose-Einstein correlations in Jet- 
set. The implementation treats them as a classical force, which violates energy-momentum 
conservation. The rescaling applied to restore the total energy and momentum, however, twists 
the event shape variables and the model description becomes worse. Multijet rates for larger 
I/cut are reduced by up to 20% and the tails of the thrust and minor distributions are decreased 
by 5-10%. The amount of particles with low rapidities is depleted by ~5%. A small but sig- 
nificant improvement is observed for small p±ont- The wave structure visible for p_Lout< 0.8 
GeV vanishes when well tuned BE parameters are used and the p_Lout distribution here can be 
perfectly described Studies on a modified implementation, which also moves unlike-sign 
pairs to avoid rescaling (additional e parameter) improves the situation but the description of 
the p° mass shift is in the wrong direction (positive). 

Another new simulation, based on the area spanned by the string, is in preparation. A first 
result with a toy Monte Carlo predicts that the reconstructed A should be 2 for 7r°, when A = 1 



is used for event generation |113 



At first glance, the experimental results are different, A ~ 1 for corrected direct measure- 
ments (DELPHI) and A ~ 2 for an extraction tuning the Jetset model. However the following 
differences must be kept in mind. The use of track mixing for a reference distribution tends 
always to give lower A values than the use of the unlike-sign meson sample. The uncorrected 
values for DELPHI are lower than for the other experiments. For kaons, corrections are es- 
timated for c- and b-decays only, but not for strong decays. ALEPH has used daughters of 
resonances wider than F = lOOMeV/c^ as prompt pions, excluding the K* which seems not to 
be affected by Bose-Einstein correlations. Ignoring this and correcting OPAL for the 7]' rate 
would bring the values down to A = 1.7 in these two analyses. 

On the model side more understanding is needed of how to include the correlation without 
twisting the event shape distribution. The new e parameter is a first step but there is no real 
success yet. Taking the decay amplitudes, i.e. string area, may be another promising approach. 
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3 Extrapolation to LEP 2 Energies 



A question of interest for LEP 2 is that of how well the characteristics of QCD events are 
understood at large energies. By QCD events, it is here meant those that are produced through 
the s-channel decay of a Z°/7* into quark and gluon jets. This question is of interest because 
W"'"W~ events lead to multi-jet states for which one of the principal backgrounds will be 
QCD events, because QCD events will also form a principal source of background for higgs, 
chargino and other particle searches, and because QCD events will be interesting in their own 
right as a means to test perturbation theory in a regime with particularly small hadronization 
uncertainties. The principal tools to test how well QCD event characteristics are understood 
are Monte Carlo generators. The main generators, Ariadne, Cojets, HERWIG and Pythia, 
have been tuned by the LEP experiments or by the Monte Carlo authors to describe global 
features of hadronic TP data. In many cases, the generators have proven able to describe detailed 
features of these data as well. It is thus relevant to extrapolate the predictions of the QCD 
generators to LEP 2 energies and to compare their level of agreement for distributions likely to 
be of importance at LEP 2. In this section, such an extrapolation and comparison is presented. 

For this study, members of each of the LEP experiments generated Monte Carlo event 
samples at Ecm= 175 GeV using parameter sets determined within their Collaboration. The 
Monte Carlo parameter sets used at LEP 1 are continually revised in order to yield as accurate 
a description of the Z° data as possible. Therefore, the parameter sets employed for this study 
do not necessarily represent official versions which will be published by the Collaborations. 
The parameter sets used for Ariadne, HERWIG and Pythia are given in tables For 
Cojets, L3 and OPAL results were made available using the parameter values given in table |^. 
There are numerous parameters and strategies involved in the optimization of the parameters. 
Comparison of the results obtained using the parameter sets of the different Collaborations 
therefore provides a systematic check of effects associated with the optimization choice. Samples 
of 100,000 events were generated without initial-state photon radiation or detector simulation, 
treating all charged and neutral particles with mean lifetimes greater than 3 • 10~^° s as stable. 

The following distributions were examined using charged particles only: 

1. charged particle multiphcity, rich, 

2. scaled particle momentum, Xp = 2p/Ecm, 

3. component of particle momentum in the event plane, p±in, and 

4. component of particle momentum out of the event plane, p±ont- 

The event plane was defined by the two vectors associated with the two largest eigenvalues of 
the Sphericity tensor. 

The following distributions were examined using both charged and neutral particles: 



1. Thrust, T [|Tl4 | 



2. Thrust major, Tmajor [|115|| 



3. Thrust minor, Tminor 115 



4. jet rates Rn defined using the k± jet finder [116|, 



23 



Parameter 


Name 


JJeiault 




UilvLr ill 


T Q 


OTO A T 


A 

Alla 


TO A TD A / 1 \ 
PAKA(ij 


U.zzU 


0.218 


0.237 


0.220 


0.200 


p_L cutoff 


TO A TD A i 0\ 

PAKA^dj 


U.oU 


0.58 


0.64 


1.00 


1 nn 

1.00 


Fragmentation function 


TV /TCT^ T / 1 1 \ 
Mbi J(ilj 


4 


3 


3 


3 


4 


Baryon model option 


TV /TCT^ T / 1 0\ 
Mbl J(l/j 


2 


o 
Z 


3 


o 
Z 


Z 


P(qq)/P(q) 


PAKJ(lj 


U.IUU 


0.100 


0.096 


n 1 nn 

0.100 


n T nn 

0.100 


Tit \ iTIt \ 

V{s)/V{u) 


1~) A T~) T /0\ 

PAKJ(2j 


U.3UU 


0.300 


0.302 


0.300 


n o nn 

0.300 


(P(us)/P(ud))/(P(s)/P(d)) 


TO A TO T /O \ 

PAKJ(dj 


U.4UU 


n /inn 

0.400 


0.650 


n /inn 

0.400 


n /I nn 

0.400 


(i/dj/-^(uaij//-^(uaoj 


TO A TO TZ/I \ 


U.UoU 


U.UoU 


U.U 


U.UoU 


U.UoU 


P(5=l)d,u 


PARJ(ll) 


0.500 


0.500 




0.500 


0.500 




PARJ(12) 


0.600 


0.600 




0.600 


0.600 


^(5=l)c,b 


PARJ(13) 


0.750 


0.750 




0.750 


0.750 


Axial, P(^=0,L=1;J=1) 


PARJ(14) 


0.000 


0.000 




0.100 


0.000 


Scalar, P(5=1,L=1; J=0) 


PARJ(15) 


0.000 


0.000 




0.100 


0.000 


Axial, ViS=l,L=l;J=l) 


PARJ(16) 


0.000 


0.000 




0.100 


0.000 


Tensor, V{S=1,L=1;J=2) 


PARJ(17) 


0.000 


0.000 




0.250 


0.000 


Extra baryon suppression 


PARJ(19) 


1.000 


1.000 


0.500 


1.000 


1.000 




PARJ(21) 


0.360 


0.354 


0.390 


0.500 


0.370 


extra r] suppression 


PARJ(25) 


1.000 


1.000 


0.650 


0.600 


1.000 


extra rj' suppression 


PARJ(26) 


0.400 


0.400 


0.230 


0.300 


0.400 


a 


PAR J (41) 


0.300 


0.500 


0.391 


0.500 


0.180 


b 


PARJ(42) 


0.580 


0.810 


0.850 


0.650 


0.340 




PAR J (54) 


-0.050 


-0.050 


-0.0378 


-0.030 






PARJ(55) 


-0.0050 


-0.0060 


-0.00255 


-0.0035 





Table 3: Optimized parameter sets for Ariadne, version 4.06 (for ALEPH, version 4.05), 
from the LEP Collaborations. The parameters listed are those which were changed from their 
default values by at least one of the groups. The Ariadne events were generated using Pythia 
version 5.7 to describe the hadronization and hadron decays. The DELPHI Collaboration 
implements its own procedure to specify the relative rate at which mesons are produced in 
different multiplets [Q, in place of the Pythia parameters PARJ(11)-PARJ(17). 
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Parameter 


Name 


Default 


ALEPH 


DELPHI 


L3 


OPAL 


Amlla 


QCDLAM 


0.180 


0.149 


0.163 


0.170 


0.160 


Cluster mass parameter 1 


CLMAX 


3.35 


3.90 


3.48 


3.20 


3.40 


Cluster mass parameter 2 


CLPOW 


2.00 


2.00 


1.49 


1.45 


1.30 


Effective gluon mass 


RMASS(13) 


0.750 


0.726 


0.650 


0.750 


0.750 


Photon virtuality cutoff 


VPCUT 


0.40 


1.00 


0.40 


0.50 


0.40 


Smearing of cluster direction 


CLSMR 


0.00 


0.56 


0.36 


0.00 


0.35 


Weight for decuplet baryons 


DECWT 


1.00 


1.00 


0.77 


1.00 


1.00 


s quark weight 


PWT(3) 


1.00 


1.00 


0.83 


1.00 


1.00 


diquark weight 


PWT(7) 


1.00 


1.00 


0.74 


1.00 


1.00 



Table 4: Optimized parameter sets for HERWIG, version 5.8, from the LEP Collaborations. 
The parameters listed are those which were changed from their default values by at least one 
of the groups. 



normalized heavy jet mass for events divided into hemispheres by the plane perpendicular 
to the Thrust axis, MheavyZ-^cm ||117|| , 



6. 
7. 
8. 

9. 
10. 
11. 

12. 



normalized difference between the heavy and light jet masses, M^^q/ 



total jet broadening, Bt ||118 



wide jet broadening, B^/ |[L18 
Sphericity, S [ p.l9| ]. 



Aplanarity, A [120 



cos6'nr| ||121|| , with four-jet 



the modified Nachtmann-Reiter four-jet angular variable, 
events defined using the k±_ jet finder with ?/cut=0.01, and 

the cosine of the angle between the two lowest energy jets in the four-jet events, cos 034. 

. were examined function of Ecm. 



In addition, the mean values of nch, Tmajor and Tj 



The results for (rich), T, Tmajor and Tminor as a function of Ecm are shown in Fig. |^ For 
those cases in which the results of at least three Collaborations are similar to each other, the 
Monte Carlo predictions are shown as shaded or hatched bands. The widths of the bands show 
the maximum deviations between the results found by the different Collaborations. The widths 
of the bands are generally much larger than the statistical uncertainties. In a few cases, the 
Monte Carlo prediction obtained by one of the Collaborations differs significantly from those 
obtained by the other three groups and is shown as a separate curve. The Cojets predictions 
are likewise shown as separate curves for purposes of clarity. The results found by the four 
LEP experiments are labelled A, D, L and O in the figure legends. 

Representative measurements from PEP, PETRA, TRISTAN and LEP 1 are included in 
Fig. ^ For -E'cm= 175 GeV, an indicative "data point" is also shown, which is taken to be 
equal to the mean of the Pythia predictions from the four groups. The size of the symbol for 
the LEP 2 point is larger than the statistical uncertainty for 10 000 QCD events. Systematic 
terms were generally found to dominate the statistical ones for the experimental measurements 
shown in Fig. 0. The total experimental uncertainties at 175 GeV can therefore be expected to 
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Parameter 


Name 


Default 


ALEPH 


DELPHI 


L3 


OPAL 


Fragmentation function 


MSTJ(ll) 


4 


3 


3 


3 


3 


Baryon model option 


MSTJ(12) 


2 


2 


3 


2 


2 


Azimuthal correlations 


MSTJ(46) 


3 





3 


3 


3 


P(qq)/P(q) 


PARJ(l) 


0.100 


0.095 


0.099 


0.100 


0.085 


P(s)/P(u) 


PAR J (2) 


0.300 


0.285 


0.308 


0.300 


0.310 


(7'(us)/P(ud))/(P(s)/P(d)) 


PARJ(3) 


0.400 


0.580 


0.650 


0.400 


0.450 


(l/3)P(udi)/P(udo) 


PAR J (4) 


0.050 


0.050 


0.070 


0.050 


0.025 


r{S=l)a,u 


PARJ(ll) 


0.500 


0.550 


— 


0.500 


0.600 




PARJ(12) 


0.600 


0.470 


— 


0.600 


0.400 


P(^=l)c,b 


PARJ(13) 


0.750 


0.600 


— 


0.750 


0.720 


Axial, V{S=0,L=1;J=1) 


PARJ(14) 


0.000 


0.096 


— 


0.100 


0.430 


Scalar, V{S=1,L=1;J=0) 


PARJ(15) 


0.000 


0.032 


— 


0.100 


0.080 


Axial, V{S=l,L=l;J=l) 


PARJ(16) 


0.000 


0.096 


— 


0.100 


0.080 


Tensor, V{S=1,L=1;J=2) 


PARJ(17) 


0.000 


0.160 


— 


0.250 


0.170 


Extra baryon suppression 


PARJ(19) 


1.000 


1.000 


0.500 


1.000 


1.000 




PARJ(21) 


0.360 


0.360 


0.408 


0.399 


0.400 


extra rj suppression 


PARJ(25) 


1.000 


1.000 


0.650 


0.600 


1.000 


extra rj' suppression 


PARJ(26) 


0.400 


0.400 


0.230 


0.300 


0.400 


a 


PAR J (41) 


0.300 


0.400 


0.417 


0.500 


0.110 


b 


PARJ(42) 


0.580 


1.030 


0.850 


0.848 


0.520 


ec 


PAR J (54) 


-0.050 


-0.050 


-0.038 


-0.030 


-0.031 


eb 


PARJ(55) 


-0.0050 


-0.0045 


-0.00284 


-0.0035 


-0.0038 


Alla 


PARJ(81) 


0.290 


0.320 


0.297 


0.306 


0.250 


Qo 


PARJ(82) 


1.000 


1.220 


1.560 


1.000 


1.900 



Table 5: Optimized parameter sets for Pythia, version 5.7, from the LEP Collaborations. 
The parameters listed are those which were changed from their default values by at least one 
of the groups. The DELPHI Collaboration implements their own procedure to specify the 
relative rate at which mesons are produced in different multiplets 0, in place of the Pythia 
parameters PARJ(11)-PARJ(17). 



Parameter 


Name 


Default 


L3 OPAL 


h 


FRAL0G(2) 


46.6 


100.0 46.6 




FRAL0G(4) 


1.52 


2.10 1.52 


K 


FRAL0Q(2) 


30.5 


43.0 30.5 




FRAL0Q(4) 


1.52 


2.10 1.52 



Table 6: Optimized parameter sets for Cojets, version 6.23, from the L3 and OPAL Collab- 
orations. The parameters listed are those which were changed from their default values by at 
least one of the groups. 
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Figure 7: The mean values of rich, Thrust T, Tmajor and Tminor predicted by Ariadne, 
Co JETS, HERWIG and Pythia as a function of Ecm in comparison with measurements from 
PEP, PETRA, TRISTAN and LEP 1. The LEP 2 point is indicative only, based on the Pythia 
prediction. The total uncertainty expected at LEP 2 assuming 10 000 QCD events is smaller 
than the symbol size. 
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be comparable to those found for the LEP 1 data. 

From the distribution of {rich) versus E^m (Fig. 0(a)), it is seen that, with the exception of 
the L3 Ariadne curve, the predictions of Ariadne, HERWIG and Pythia are similar. The 
widths of the Ariadne, HERWIG and Pythia bands are narrow for energies at and below the 
Z° mass, showing that the results from the four Collaborations are in close agreement (with the 
exception of the L3 curve for Ariadne). For energies above about 150 GeV, the HERWIG band 
becomes broader, indicating that there is some divergence in the predictions obtained by the 
different groups. From Fig. |^(a) it is also seen that COJETS predicts a substantially larger value 
of {rich) than the other models for energies above the Z° mass. This difference is suggestive of 
coherence effects in the parton shower, which are absent in CoJETS but present in the other 
three models. Coherence reduces the mean soft gluon multiplicity in the parton shower. It is 
generally expected that coherence will lead to a reduction in the mean hadron multiplicity as 
well. Thus, a measurement of {rich) at LEP 2 could help to establish the existence of coherence 
phenomena in the data. 

Figs. |^(b)-(d) show the corresponding distributions for the T, Tmajor and Tminor variables. 
Again, Ariadne, HERWIG and Pythia are seen to exhibit similar behavior. Cojets agrees 
well with the other models for T, but lies below them for Tmajor and above them for Tmmor in 
the LEP 2 energy range. Thus the jets from CoJETS are less oblate than those from Ariadne, 
HERWIG or Pythia. (The Oblateness O of an event is given by O = Tmajor — ^minor-) The 
differences between Cojets and the other three models become larger as Ecm increases. 

In Fig. H, the Monte Carlo predictions for rich, Xp, p±in and p±ont at 175 GeV are shown. 
The corresponding results for T, Tmajor, Tminor and for Mheavy/-£^cm, Miifr/Ecm, Bt and Bw, 
and for S, A, |cos6'J^j^| and cosa34 are shown in Figs. H, 0, and [IT], respectively. Overall, 
the models are seen to be in general agreement with each other. Some of the more notable 
exceptions to this agreement are discussed below. 

1. A striking difference is observed between Co jets and the other models for the rich and 
P±m distributions (Figs. ||(a) and (c)). Smaller but visible differences are observed between 
Cojets and the other models for a number of the other distributions as well. At the Z" 
mass, these differences between Cojets and the other models are either not present or 
are much smaller. This implies that the energy scaling behavior of Cojets differs from 
that of Ariadne, HERWIG and Pythia. 

2. For HERWIG, the Xp distribution is much harder using the L3 parameter set than it is 
using the parameter sets of the other Collaborations (Fig. Kb)). This feature is also 
observed at the Z° energy. The primary reason for this difference between L3 and the 
other groups is the different treatment of the parameter CLSMR (see table ^). 

3. From Fig. |^(d), it is seen that the three-jet rate from Pythia is significantly larger than 
that of the other models for ?/cut values below about 0.02. Correspondingly, the two jet 
rate from Pythia is smaller. This difference is also observed at Ecm= 91 GeV. From this 
same figure, Cojets is seen to predict a three-jet rate which is smaller than that of the 
other models: this last difference is not observed at LEP 1 energies. 

4. Cojets exhibits a clear deviation with respect to the predictions of the other models for 
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Figure 8: Comparison of the predictions of QCD event generators at Ecm= 175 GeV. 
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Figure 9: Comparison of the predictions of QCD event generators at Ecm= 175 GeV. 
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the jet mass distributions, MheavyZ-Ecm and MdifrZ-Ecm (Figs, ^{a.) and (b)). Less of a 
deviation is present for the jet broadening variables, Bt and Bw (Figs. |lO|(c) and (d)). 
This suggests that these last two variables may be less subject to uncertainties related to 
the modelling of QCD and hadronization than the first two variables. 
5. CoJETS and HERWIG are seen to exhibit a somewhat flatter distribution in |cos6'^f^| than 

Ariadne and Pythia (Fig. |ll](c)). 
The general conclusion that can be drawn from this study is that there is relatively little 
uncertainty in the predictions of QCD generators for event characteristics at LEP 2. Such basic 
features of events as charged multiplicity. Thrust and Oblateness are described in an almost 
identical manner by Ariadne, HERWIG and Pythia. Only CoJETS deviates significantly 
from the predictions of the other models. On the other hand, there is modest disagreement 
between the models for variables which require use of a jet finding algorithm: (Fig. §(d)) 
and |cos6'^j^| (Fig. 0(c)). This could have some implication for the W mass determination 
based on the reconstruction of jets. 
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Figure 10: Comparison of the predictions of QCD event generators at Ecm= 175 GeV. 



32 



Figure 11: Comparison of the predictions of QCD event generators at Ecm= 175 GeV. 
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4 Monte Carlo descriptions 



In this section we have collected brief descriptions for the main QCD generators and other pieces 
of QCD code. These writeups are intended to introduce the main physics ideas and give further 
references to manuals and codes — a full coverage of all physics and programming aspects is 
excluded for space reasons. The compilation below should be rather complete for programs 
intended for the main QCD-related processes, such as 7*/Z° and W pair production. Special 
emphasis is put on HERWIG and P YTHIA / Jetset, which have been used extensively at LEP 1 
and are equipped with a simulation both of electroweak and QCD aspects. A few programs 
include QCD aspects but have still been judged to better belong elsewhere, e.g. Phojet is a 
77 physics generator and Isajet is mainly of interest (in the e"'"e~ sector) as a supersymmetry 
generator. 



4.1 ARIADNE 



Basic Facts 



Program name: 

Version: 

Author: 



Ariadne 13 



Program size: 
Program location: 



4.08 of 30 November 1995 
Leif Lonnblad 

NORDITA, Blegdamsvej 17, 
DK 2100 Copenhagen 0, Denmark 
Phone: + 45 - 35325285 
E-mail: leif@nordita.dk 
12853 hues 

http:/ /suryall.cern.ch/users/lonnblad/ariadne/ 



The Ariadne program implements the Dipole Cascade Model (DCM) for QCD cascades 
||122|| . In this model the emission of a gluon g^^ from a qq pair created in an e"'"e~ annihilation 
event can be described as radiation from the colour dipole between the q and q. A subsequent 
emission of a softer gluon gg can be described as radiation from two independent colour dipoles, 
one between the q and g^ and one between and q, neglecting the contribution from the qq 
dipole, which is suppressed by 1/Nq. Further gluon emissions are given by three independent 
dipoles, etc. In this way, the end result is a chain of dipoles, where one dipole connects two 
partons, and a gluon connects two dipoles. This is in close correspondence with the Lund string 
picture, where gluons act as kinks on a string-like field stretched between the qq pair. 

This formulation of the partonic cascade in terms of colour dipoles means that the coherence 
effects, handled by introducing angular ordering in conventional parton showers, is correctly 
taken into account. Also, the DCM has the advantage that the first gluon emission is done 
according to the correct first-order matrix element, so that an explicit matching procedure like 
the ones introduced in HERWIG and Jetset is not needed. 
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Although the model has been developed a lot since the last LEP workshop, much of this 
development has been related to the description of Deep Inelastic Scattering and hadron-hadron 
collisions and will be described in some detail in the report from the 77 generator working group. 
Here only aspects relevant to e"'"e~ annihilation will be discussed. 

The basic DCM only describes gluon emission, and the process of splitting a gluon into 
a qq pair has therefore been added according to [|123|| . Although this procedure reproduces 



fairly well the amount of secondary cc production observed at LEP ||124|| , there has been some 



criticism [|125|| that the model may be overestimating the phase space available for this process. 



Therefore an extra restriction of this phase space suggested in [ p.25|| has been implemented as 
an option in the last versions. 

The radiation of photons from quarks is handled by allowing the process of emitting a photon 
from the electro-magnetic dipole between the original qq to compete with the gluon emission 
from the colour dipoles ||126|| . This competition is governed by the ordering in transverse 
momenta of the emitted gluons/photons, which is different from Jetset and HERWIG, where 
virtuality and angle, respectively, is used for ordering. 

In the latest version, a scheme for colour reconnections has been added to the program. The 



model is described fully in [|127|] and briefly in section |4.7| . Unfortunately, the manual included 
in the code distribution has not yet been updated to describe this new feature, and users who 
want to try it are advised to contact the author by e-mail before doing so. 

Since Ariadne only handles the perturbative QCD cascade in an event, it has to be inter- 
faced to the Pythia/ Jetset programs for generation of the hard sub-process, the hadroniza- 
tion and the particle decays. Such an interface is included in the code, and only very minor 
changes to the steering program is needed to replace the parton showers in Pythia/ Jetset 
with the dipole shower in Ariadne for any type of process. In a typical steering program for 
running Pythia, the changes needed are as follows. 

• Immediately before the call to PYINIT there should be inserted a call to ARTUNE( '4.07' ) 
to set up default parameters in Ariadne and Pythia/ Jetset, followed by a call to 
ARINIK 'PYTHIA') to initialize the Pythia interface. To change the default behavior 
of Ariadne, changes may be made to the ARDATl common block between the calls to 
ARTUNE and ARINIT. 

• Immediately after a call to PYEVNT, a call to AREXEC should be made to perform the actual 
dipole cascade. If Pythia is set up to handle fragmentation and decays in the PYEVNT 
call, this is now handled in AREXEC instead. 

Sample programs for how to do this is included in the distribution. The distribution also 
includes a subroutine AR4FRM which is an interface to four-fermion generators according to 
the standard presented in section ^.3| . Except for what is needed to run Pythia/ Jetset, no 
additional software is required to run Ariadne. 
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4.2 COJETS 



Basic Facts 
Program name: 
Version: 
Author: 



COJETS ME 



Program size: 
Program location: 



CoJETS 6.23 of 10 February 1992 

Roberto Odorico 

Department of Physics 

University of Bologna 

Via Irnerio 46, 1-40126 Bologna, Italy 

Phone: + 39 - 51 - 24 20 18 

E-mail: odorico@bo.infn.it 

19742 lines 

http: / / www, bo.infn.it / preprint / odorico.html 



CojETS simulates electron-positron annihilation into jets of hadrons. (It also simulates 
proton-proton and antiproton-proton interactions.) The simulation is based, at the parton 
level, on the standard model with perturbative QCD treated in the leading-log approximation. 
QED radiation off beam particles is treated according to the BKJ program. Partons from parton 
showers are independently fragmented into jets of hadrons according to a Field-Feynman model 
extended to include heavy quarks and baryons and modified in the generation of soft particles. 
Gluons are fragmented as a pair of light quark and antiquark jets of opposite random flavors, 
each one having half the energy of the gluon and its same direction and with fragmentation 
parameters distinct from those of quark jets. Jet non-perturbative masses are limited by bounds 
originated by an approximate treatment of global phase-space effects at the multi- jet level. 
A previous version of the program, CoJETS 6.12, in which quarks and gluons share the same 
fragmentation model is also available. The jet fragmentation model adopted goes hand in hand 
with the setting of the minimum parton-mass cutoff to a value of 3 GeV, which is substantially 
larger than those used in string- and cluster-based fragmentation models. 

The output common block, containing the generated particle stream, has the standard 
/HEPEVT/ format, with PDG codes used for particles. 

CoJETS is maintained with the PATCHY code management system. The appropriate FOR- 
TRAN77 codes are obtained by means of suitable pilot patches. The program file also includes 
the documentation. 

Recently, the program has been mainly used to check the relevance of evidence for string 
fragmentation, parton coherence and quark/gluon differences in jet fragmentation. Its usage to 
study the signal/background enrichment for bottom non-leptonic decay events by means of neu- 
ral networks has shown that differences in internal correlations between signal and background 



are fuzzier in Cojets than in Jetset 



Thus for LEP 2 Cojets could be useful when 



studying ways of disengaging events with W pairs decaying non-leptonically from background. 

Cojets had its fragmentation parameters sensibly tuned to reproduce basic experimental 
distributions. So far the tuning has been done by the author, mimicking experimental apparatus 
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effects but without a proper GEANT simulation. That can be done by the user by means of 
the program TUNEMC, based on Minuit's Simplex algorithm (a more advanced version of the 
program is in preparation). 

Programs Co JETS 6.12 and TUNEMC can be found at the same WWW URL as Co jets 
6.23. 



4.3 HERWIG 



Basic Facts 

Program name: 

Version: 

Authors: 



Program size: 
Program location: 



HERWIG 14 



HERWIG 5.9 from 1 January 1996 

G. Marchesini^ B.R. Webber^, G. Abbiendi^ I.G. Knowles^ 
M.H. Seymour^, L. Stanco^ 

1, Dipartimento di Fiscia, Universita di Milano. 

2, Cavendish Laboratory, University of Cambridge. 

3, Dipartmento di Fisica, Universita di Padova. 

4, Department of Physics and Astronomy, University of Glasgow. 

5, Theory Division, CERN. 

E-mail: webber@hep.phy.cam.ac.uk, knowles@v6.ph.gla.ac.uk, 
Seymour @surya 1 1 . cern . ch . 

15500 hues ^ 

http: / / suryal 1 .cern.ch / users/ seymour/herwig/ 



4.3.1 Introduction 

HERWIG (Hadron Emission Reactions With Interfering Gluons) is a large, multipurpose Monte 
Carlo event generator which has been extensively used at LEP 1. Version 3.2, as applied to 



e"'"e annihilation, was described in detail for the LEP 1 workshop ||130||, here we concentrate 



principally on program developments and new aspects of relevance to LEP 2 physics. 

QCD Monte Carlo event generators utilize the fact that any hard scattering processes can 
be factorized into separate components at leading twist. These are: the hard sub-process 
itself; perturbative initial and final state showers; non-perturbative hadronization; resonance 
decays; and beam remnant fragmentation. In HERWIG great emphasis is placed on making 
available a very sophisticated, partonic treatment of the calculable QCD showers. In contrast 
the description of the at present uncalculable hadronization and beam remnant components is 
in terms of very simple models. Since HERWIG contains many hard sub-processes and supports 
all combinations of hadron, lepton and photon beams this allows the physics of many types 
of particle collisions to be simulated in the same package. In view of the universality of the 
factorized components that build up HERWIG events this allows experience gained at HERA 
and the TEVATRON, for example, to be made directly available to LEP physicists. 
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Since version 3.2 was released the HERWIG code was reorganized to isolate the shower, 
cluster hadronization and unstable particle decay routines. This modularity facilitates the 
creation of hybrid programs in which sections of code are replaced with interfaces to other 
Monte Carlo programs. To identify HERWIG code the names of all options statements now begin 
HW****. The /HEPEVT/ standard proposed in [|131|| is also now used throughout the program in 
DOUBLE PRECISION. The random number generator has been upgraded to a I'Ecuyer's algorithm 
as recommended in [[L32|| . Discussion of physics changes are contained in the following sections. 



4.3.2 Hard Sub-processes 

An extensive range of hard sub-processes are available in the HERWIG program allowing the 
full spectrum of standard model LEP 2 physics to be simulated. These are illustrated in table 
for a complete listing see the program release notes in the text file HERWIGnm.DOC. 

The matrix elements used for the continuum processes IPR0C=100-153 now allow for arbi- 
trary polarization of the lepton beams, an additional Z' including complete 7*/Z° interference 
and full mass effects (IQ 7^ 0). When the ZPRIME= . TRUE . option is set the Z' weak couplings 
used are taken from the arrays AFCH(*,2) and VFCH(*,2), see HWIGIN for details. The arrays 
Q/V/AFCH(*, 1) are used consistently throughout the program for the standard model electric, 
weak vector and axial- vector fermion couplings. A running electromagnetic coupling aem(Q^) 
is used for internal photons |p.33|| with the hadronic part taken from [|134|| . It is normalized to 



the Thomson limit [Q^ = 0) value ALPHEM. Process 107 is included to facilitate q/g studies; 
in analogy with quarks the gluons are given a 1 -|- cos^ 6 distribution. The difference between 
processes 100-106 and the original 120-126 lies in the treatment of hard gluon emission. The 
massless matrix element matching scheme used by HERWIG is discussed under parton showers. 

Specialist Matrix Element programs for W'''W~/Z°Z° production, more properly four fermi- 
on generators, employ the full set of gauge invariant diagrams. In comparison HERWIG only 
includes the subset of diagrams containing W+W~ ("CC03") or Z°Z° ("NC02") pairs but does 
provide realistic hadronic final states. These matrix elements are taken from the program of 
Kunszt ||135|| and correctly include spin correlations in the gauge boson decays. Additionally a 
model for colour re-arrangement within the context of HERWIG is available, see [4.7| for details. 
The decays of the vector bosons are controlled via the array MODBOS, as detailed in HWIGIN, and 



include spin correlations. Please also see the detailed prescription, discussed in section |5]^, for 
interfacing HERWIG to specialist four fermion Monte Carlos. 

At LEP 2 the principal Higgs production mechanism is the Bjorken process, IPR0C=300+ID, 
where one or both Z°'s may be off-shell, Z"*-*-* Z°'-*''h°; also available is vector boson fusion, 
IPR0C=400-|-ID. In both cases the exact leading order matrix element is used in the improved 



s-channel approximation ||136|| . At LEP 2 a discoverable Higgs would be narrow; in the program 
the actual mass used is taken from within the range [Mh — GAMMAXT^rPh, Mh + GAMMAX^irrh] (default 
GAMMAX=10) using a 'Breit-Wigner' distribution corrected for an energy dependent width. The 
event weight is the product of the cross-section (in nb) multiplied by the branching ratio to the 
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channel specified by ID. The Higgs' partial widths are calculated in HWDHIG: the quark decay 
channels include next-to-leading logarithmic corrections and the vector boson decay modes 
allow for off-shell WW/ZZ pairs. 

The cross-sections for 77 interactions rise with cm. energy to become the commonest 
physics processes at LEP 2. When considering hadronic final states each photon may be viewed 
as interacting either as a point-like particle or as being resolved into constituent (anti-) quarks 
and gluons. This leads to three basic sets of hard sub-processes (zero, singly or doubly resolved), 
a division adopted in the wide selection of sub-processes made available in HERWIG. Note that 
this separation is in fact artificial and all three components must be combined to obtain the full 
cross-section. Discussion of these processes can be found in the HERWIG description provided 
in the gamma-gamma section of this report. 



4.3.3 Initial State Radiation 



In e^e~ scattering real photons are radiated from the incoming lepton lines. At LEP 1 any 
effects were mitigated against by the penalty involved in going off the Z° resonance. However 
photon bremsstrahlung is expected to be an important feature at LEP 2 energies where the 
basic cross-section typically rises as s decreases. HERWIG uses an electron structure function 
approach to write the total cross-section for a process as: 



fl rl 

dxi / dx2f^{xi)f^{x2)a{xiX2s) (2) 
) Jo 

Employing a natural choice of variables, r = xiX2{= s/s) and x = xi, this can be written: 

/;dr*M/;!^/;W/;(I) (3) 



where r > T (TMNISR) is a physical cut-off used to avoid the 1/s pole in the cross-section's 
photon exchange term. 

The actual structure function used is the second order solution to the full Altarelli-Parisi 
equation with exponentiated coefficients: 



f:{x) = /5(l-^r'exp{/3|(l + |)}(l-/3^£ 



P 



71 



log (^^ 
[ml , 



1 



+ 



;i + a;)[(l + x)2 + 31oga;] 



4 logx 

1 — X 



(4) 



This is equivalent^ to the expression, eqs. (58,60), given on p. 34 of This means that 
the single photon emission allowed for in HERWIG gives equivalent energy and p-^ spectra as 

discrepancy in the coeflicient of the tt^ term, a factor 2 too large, is beheved to be their typographic 

error. 
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IPROC 



Process 



^* /Z^ /Z' Continuum Processes 



100+IQ 
107 

110+IQ 
120+IQ 
127 

150+IL 



e ' e 

c+e" 

e"'"e" 

e+e" 

e''"e" 

e+e" 



qq(g) 

gg(g) 
qqg 
qq 
gg 

u 



IQ=l-6: q = d, . . . , t; IQ=0: all flavours 

IQ as above, includes masses exactly 
IQ as above T no correction to hard 
J gluon branching 

1L=2,3: t = 



Di-Boson Production 



W^/Z*^ decays controlled 
by MODBOS 



200 
250 



e ' e 
e+e" 



} 



+h^ 

+hO 

+h° 

+h° 

+h° 



Bjorken process: e^'^e — > Z h 



300+IQ 

300+lL 
310,311 
312 
399 



qq 

u 

w+w-,z°z° 

77 

anything 



IQ as above 
IL=1,2,3: e^e,fx,T 



Vector Boson Fusion 



e+e-h^ 



400+ID 



e ' e 



ID as 1PROC=300+1D 



Zero Resolved Gamma- Gamma: c^c 



qq/££/W+W- 



(c+e )77 



500+lD 



77 



1D=0-10 as for 1PROC=300+1D 



Gamma-H^ Fusion: e^e 



-//e)7W- 



550+ID 



7W- 



qq'/ 



ID=0-9 as for IPROC=300+ID 



Doubly Resolved Gamma-Gamma 



1500 
1700+lQ 
1800 
2200 



gq 

qq 



g&qg ^ 
QQ, gQ 
7q, gg - 

77: gg - 



qg, etc. 

gQ, etc. 
* 7g, etc. 

^ 77 



31 0{a'^) two-to-two QCD scatterings 

16 0(as) heavy quark production processes 

17 0{asacm, «s<^cm) direct photon processes 

„2 



3 Oia. 



cm' s cm/ 



di-photon processes 



Singh' Resolvcxl G annua- G annua 



5000 
5100+lQ 
5200+lQ 
5500 

5510,5520 



7q gq, 7g qq, etc. 3 ©(agCKem) dijet processes 
7g QQ heavy flavour pair production, IQ as above 

7Q ~^ gQ, qQ ~^ qQ heavy flavour excitation, IQ as above 

gg — > Vg, gq — s> Vq', etc. 6 0{agaem) light (u,d,s) L = meson production 

J{— S) — 0,1 mesons only 



8000 



Minimum bias soft coUision 



Gharged lepton Deep Inelastic Scattering 



9000+lQ 
9010+lQ 



eq eq, cq eq 
eq z/eq', eq t^eq' 



NG DIS on flavour IQ as above 
GG DIS on flavour IQ as above 



10000+IP 



As IPROG=IP but with suppressed SUE 



Table 7: The principal HERWIG hard sub-process of importance for LEP 2 physics. In QGD 
scatterings IHPRO labels the actual sub-process, allowing for colour decomposition, generated. 
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multiple photon emission. In the soft photon limit fc{x) simplifies significantly to the following 
form, used to efficiently generate the {xi} via importance sampling: 

f!{x,Q')^Pil-xr-' (5) 

In practical situations one has: acm/7r <^ /? -C 1, so that f^ix) has an integrable singularity in 
the soft photon limit, x ^ 1; 1—x is the energy fraction carried by the photon. To regularize this 
divergence HERWIG employs a resolution parameter, x < X, called ZMXISR (default 1 — 10~^), 
and includes a fraction of events with no emission, so that: 

f:ix) ^ fix) = e(X - x)fix) + 6(1 - x)(l - Xf (6) 

Here X is only an internal parameter and unphysical in the sense that cross-sections should not 
depend on it. Observe that even for 1 — X = O{10~^) the non-emission probability is ~ 45% 
at LEP 2. Note also setting ZMXISR=0 has the effect of switching off the initial state photon 
radiation. 

After the emission of a photon the electron entering the hard sub-process is off-shell. In 
HERWIG its negative virtuality is selected from a logarithmic distribution, dg^/g^, bounded in 
magnitude by eq.(H). 

Allowing for the virtualities of the electron lines and treating a: as a lightcone momentum 
fraction s is reconstructed as: 

s = TS-q^-q2 + 2p p (7) 

s — 

Since a is a rapidly varying function of s near the Z° HERWIG slightly shifts the {xj fractions 
to preserve s = ts. Specifically the highest p^ photon is taken to be emitted first and its 
X x' = X + (f' I s {x' is the energy fraction) so that s would be preserved in the absence of 
emission from the other lepton. The x of the lower p^ photon is then shifted so as to give 
exactly s. For simplicity the program requires photon emission to be in the forward hemisphere 
which imposes the condition: 

q'<^{l-x)s (8) 

This inequality is applied to both leptons. Note that this still allows the possibility for on 
resonance 7Z° states to be produced but only to the accuracy of the leading logarithm approx- 
imation. 

The use of the Equivalent Photon Approximation for the case of virtual photon emission in 
which it is the photon which enters the hard sub-process is again discussed in the report of the 
gamma- gamma working group. 



4.3.4 Parton Showers 



HERWIG employs highly developed parton shower algorithms to provide an accurate description 
of the perturbative QCD jet evolution. Coherence, due to leading infrared singularities [^, is 

41 



automatically included through the choice of evolution variables, ordering in which naturally 
restricts the branching phase space to an angular ordered region. Further angular screening 
due to heavy quark masses, the dead cone, is also fully included ||137[ . At each branching the 
azimuthal angles are distributed according to the eikonal dipole distribution for soft gluons ||138|| , 
including mass effects, and to the full coUinear leading logarithm accuracy for hard emission 
At large momentum fractions the coherent algorithm used also correctly describes next- 

By using a two-loop expression for as this allows the Monte 
1 



139 



to- leading contributions ||140 | 



Carlo A to be related to Aj^ 



A 



MC 



exp 



C^(67- 37r2) - lOnj 
6(llC^-2n/) 



A — » 



.5l.569A(| 



(9) 



Since the time of the LEP 1 workshop significant progress has been made in the study of final 
states involving photons ||103|| , leading to the implementation of final state photon radiation in 



HERWIG ||141|| . The momentum sharing in q ^ q7 branchings and relative rate compared to 



q — s> qg branchings are controlled by the following splitting function and Sudakov form factor: 



™^ ^ ^27r 1-z 

2 ^cm 



logA,(g^g^) = -ei^[{\og{Q/Q,)~3/Ay~{\og{Qo/Q,)-3/4y] (10) 



TT 



where, since the photon is in the final state, a fixed a^m is used (allowing analytic integration 
of the Sudakov form factor) and Qo = Qq + Q-y with Qq and the cut-offs on the quark and 
photon scales respectively. The branching 7 — > qq is expected to be small and is not included. 
Competition between the two types of quark branching is handled in the standard way. That is 
the scales at which the two types of branching attain a preselected probability of occurring 
are found, the larger is taken to occur first and if its is above (Qq + Q^y)^ it is accepted. The 
scale of any branching is bounded above by that of the last emission, irrespective of type. 
However the opening angle is bounded from above by the opening angle of the last emission of 
the same type; this is exact in the case when azimuthal photon-gluon correlations are integrated 
out. 

Due to the choice of evolution variables in HERWIG, configurations in which a very hard 
gluon or photon recoils against the qq pair are not generated by the showering algorithm, that is 
a 'dead zone' exists [|141|| . This is particularly important in the photon case due to the relative 



ease with which they can be identified in the final state. The HERWIG solution is to find what 
fraction of events are missing by integrating the three parton matrix element over the dead 
zone and then add back this fraction starting the evolution from a correctly distributed qqg/ 7 



configuration. The algorithm of ||142|| is used to exactly include initial/final state correlations 



starting from a massless qq configuration. 

The matching of a hard gluon or photon to the exact matrix element is controlled by the 
logical HARDME (default .TRUE.). Additionally there is a 'soft' matrix element correction, where 
soft here means inside the phase space region accessible to the branching algorithm; SOFTME 
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(default .TRUE.) controls the matching of the hardest emission, not necessarily the first, to the 
exact matrix element 11143 1 . 



4.3.5 Hadronization 



The basic preconfinement inspired |73[ cluster hadronization model used in HERWIG remains lit- 



tle changed from its original formulation ||144|| . The principle criterion for selecting the flavours 
and spins of the primary hadrons in the cluster two body decays is the phase space available; 
though weights PWT (1-6) , VECWT , TENWT and DECWT can be used to alter the flavour/spin com- 
positions. The cluster decays are isotropic, in their own rest frame, except when a perturbative 
quark is involved, that is one from the hard sub-process or a g — > qq splitting. If (CLDIR=1), the 
default, then the hadron containing this quark is aligned with the quark direction in the cluster 
rest frame. The main effect is to stiffen the spectrum of heavy charm and bottom hadrons. It 
is possible to partially decorrelate this direction retention using the parameter CLSMR (default 
0), the width of an exponential distribution in 1 — cos6'qh; thus increasing CLSMR increases the 
smearing. 

New parameters have been introduced to control the treatment of clusters with anomalous 
masses. CLPOW (default 2) influences the decision on whether a heavy cluster q^^qg should first 
be split in two prior to hadronization according to if its mass satisfies the inequality: 

MClpow > c-LMAX^^^o^ + (mq^ + m^J^^^^^ (11) 

Using smaller values of CLPOW leads to an increased yield of heavy clusters containing heavy 
quarks and thence to more heavy baryons; light quark clusters are affected less. The parameter 
BILIM (default 0) can be used to increase the number of relatively light bottom clusters that 
undergo a one-body decay. If Mthr is the threshold for two-body decay then the probability of 
a one-body decay becomes: 

r 1 M,i < Mthr 

'P={ ^-mtkMlL ^^^r < M,, < {l + BlLIM)Mthr (12) 

[ {1 + BlLIM)Mthr < Md 

For light quark clusters the one-body decay criterion remains equivalent to the above with 
B1LIM=0. In practice CLPOW proves more effective in controlling the spectrum of both bottom 
and charm hadrons. When one-body decays do occur a Lorentz covariant treatment is now 
used to effect the necessary momentum rearrangement. 

In the default version of HERWIG the quark-antiquark pairs which form the colour singlet 
clusters are taken to be nearest neighbour pairs, in a sense defined by the shower. However a 
colour reconnection model is now available. It is based upon minimizing the spatial sizes of 
pairs of clusters as determined from the semi-classical positions of the partons at the end of 
the showers. This model is discussed more fully in I 



More recently the number of hadrons supported has been enlarged to incorporate all L = 0, 1 
mesons (including the 0"*"*^+^ and l+(+) states) composed of d,u,s,c,b quarks and all J = 1/2 
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('octet') and J = 3/2 ('decuplet') baryons composed of the form qiq2q3 or Qqiq2, (Q = c,b). 
Should the user wish to add any new particles it is sufficient to simply specify their properties: 
name, PDG code number, mass, spin and flavour compositions in the arrays RNAME, IDPDG, 
RMASS , RSPIN and IFLAV and they will be included automatically in cluster decays. Using the 
array VTOCDK it is also possible to veto a particular hadron's production in cluster decays. 



4.3.6 Decay Tables 

The HERWIG decay routines have been largely re-written to make them more user friendly 
and to adopt the proposals made in section pl2| . Up to five body decays are supported with 
a number of standard matrix elements made available. Specific hadronic decay channels for 
B hadrons can now be included. This is in addition to the original partonic model based 
on spectator decays | |137| |; note this may involve some double counting. The production of a 



selected particle via unstable particle decays can be vetoed by specifying it in the array VTORDK; 
any branching ratio sums affected because of excluded channels are automatically reset to unity. 
The subroutine HWIODK has been added to allow the HERWIG decay tables to be inputted and 
outputted in the proposed standard format. When read in the program checks that the decay 
is kinematically allowed and does not violate electric charge conservation; if necessary the sum 
of branching ratios is reset to one. The use of this subroutine makes it simple for the users to 
adapt the provided tables for their own use. The subroutine HWMODK allows individual channels 
in the decay tables to be added or modified between events. The actual default decay tables 
themselves have also been updated to include modes at the one per mille level. 



Interfaces to the EURODEC ||145[ and CLEO p46 | B hadron decay packages are also built into 



HERWIG. The selection is made by setting BDECAY= ' EURO ' , 'CLEO' or 'HERW (the default is 
of course 'HERW'). 

The production vertices of hadrons are now calculated by HERWIG and stored using the 
VHEP array of /HEPEVT/. This is based on the particle lifetimes in the RLTIM array. A particle 
is set unstable if its lifetime is less than PLTCUT however when MAXDKL= . TRUE . all decays are 
tested in the routine HWDXLM and required to occur within a volume specified by lOPDKL else 
left undecayed. If BOMIX=.TRUE. then neutral B^^ mesons are allowed to mix before decaying. 



4.3.7 Source Code 



In addition to the WWW site quoted above copies of the HERWIG source code and supporting 
files are maintained in the following VAX directories: 

CBHEP: :DISK$THEORY: [THEORY . HERWIG] HERWIGnm . * 
FNALV: :USR$R00T2: [BWEBBER . HERWIG] HERWIGnm. * 
VXCERN: :DISK$CR: [WEBBER. HERWIG] HERWIGnm. * 
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ThefilessuppliedareHERWIGnm.COM, *.DOC, *.FOR, *.INC, *.MSG, * . SUD and * . TST. The 
command file HERWIGnm.COM runs a test job * .TST containing the main program. This uses the 
source code subroutines found in * . FOR with the declarations and common blocks in * . INC and 
default Sudakov form factors in * . SUD. Release notes are found in * . MSG and more complete 
documentation in * . DOC. 



4.4 NLLjet 



Basic Facts 

Program names: 

Versions: 

Author: 



Program size: 



NLLjet 



NLLjet 3.0 of September 1992 
Kiyoshi Kato 
Kogakuin University 

Nishi-Shinjuku 1-24, Shinjuku, Tokyo 160, Japan 

Phone: + 81 - 3 - 3342 - 1211 

E-mail: kato@sin.cc.kogakuin.ac.jp 

Tomo Munehisa 

Yamanashi University 

Takeda 4-3, Kofu 400, Japan 

Phone: + 81 - 552 - 20 - 8584 

E-mail: munehisa@top . esb . yamanashi. ac . j p 

7742 lines 



Program location: ftp.kek.jp : kek/minami/nlljet 



NLLjet is a Monte Carlo code for the generation of jet events in e"'"e~ annihilation based 
on the parton shower method. The events are parton final states in the form of a list with 
particle codes and four-momenta. Connection to the hadronization is open for the user, and a 
standard interface to Lund hadronization is provided. 

Generation of QCD jets by the parton-shower method was born of Konishi, Ukawa and 
Veneziano in 1979 as the "jet calculus" in which the method to make systematic summation 
of the collinear singularity in QCD was given. Here, the factorization of the mass singularity 
works well and the choice of physical gauge leads to a suppression of interference terms, so that 
a stochastic treatment for jets becomes possible. 

Soon after that, models of the QCD parton shower in the leading-logarithmic (LL) approx- 
imation were developed. These models are good for the description of jets in high energy. 
However, they have no chance to determine the fundamental parameter of QCD, Q;s(/i^) (or 
QCD A), because starting from any renormalization scheme in QCD, you obtain the same 
formula for physical quantities in LL approximation. This limits the analysis for the determi- 
nation of the strong coupling constant in jet phenomena only to the calculation based on the 
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QCD matrix elements. However, the Monte-Carlo simulation of jets by matrix elements is not 
appropriate for the global description of jets since it has an avoidable defect, the discontinuity 
between n- and (n + l)-parton states. 

The idea of NLLjet was spawned from observation above. In this parton-shower model, 
the collinear singularity of QCD is summed up to the next-to-leading logarithmic(NLL) order. 
All components in NLL order are computed in the MS scheme, and they are implemented in 
the model. Thus NLLjet has the potential to determine the QCD Ky^ through a comparison 
of generated events with experiments [|148|| . The basic ingredients of NLLjet are as follows: 



• Sudakov factor which is defined by the integral of the P function up to 0{a1). 

• Two-body branching by the two-body vertex function up to Oi^al). 

• Three-body branching by the three-body vertex function in 0{a1). 

• Hard cross section of the primary qqg process up to 0{al). 

• Kinematical conditions and correction terms. 

The effect of soft-gluon contribution is an important issue in perturbative QCD. In NLLjet, 
the strong coupling constant in the Sudakov factor is defined to be as{x{l — x)Q'^), and it 
corresponds to the inclusion of soft gluon resummation. The angular ordering is not introduced 
to all branchings but only to those in which the angular ordering is really required. 

The important point of the formulation beyond LL order is that each kinematical modifi- 
cation is always controlled properly through the introduction of a correction term in the NLL 
order functions. The three-body vertex functions become positive with the correction for the 
angular ordering inq— *q + g + g and g — * g + g + g and that for the momentum conservation. 
The double cascade scheme, which is necessary to recover the symmetry between q and q, also 
gives another correction term. 

The parton shower method still has a few ambiguous points which are hard to determine 
from the theoretical view point in perturbative QCD. For example, the virtuality of partons in 
final states should be less than a cutoff value, Q^. Normally, one sets it equal to Q^. However, 
sometimes better agreement with experiments is found by taking it to be 0. In this version, 
this modification is included by setting KINEM -1 parameter. 

The effect of a quark mass is only counted kinematically by replacing by + m^. 
Neither azimuthal correlations nor the parton polarization are considered. 

Essential input parameters of NLLjet are W, A, Ql, 6, and C. Here W stands for the 
center-of-mass energy. Physics should not depend strongly on the cutoff Ql and its dependence 
is to be counted as a systematic error of the theory. Parameter 6 is specific to NLL parton 
shower and it is absent in LL order. The distribution is expected to be independent of 6. 
However, detailed study shows that there is small bend at the region connected by 6. If one 
sets 6 large (~ 0.5), the events are free from the bend at the expense of the exclusion of qqg 
primary vertex. The scheme parameter C is available in order to replace /x^ in as{fi'^) from 
/i^ = to /i^ = CQ^. However, it is not possible to change C in large. In the matrix element, 
QCD is studied at = W"^ while in parton shower, it is done for = W"^ ~ Qo- 



46 



4.5 PYTHIA/JETSET 



Basic Facts 



Program names: 
Versions: 

Author: 



Program size: 
Program location: 
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Pythia and Jetset 
Pythia 5.720 of 29 November 1995 
Jetset 7.408 of 23 August 1995 
Torbjorn Sjostrand 
Department of Theoretical Physics 
University of Lund 

Solvegatan 14A, S-223 62 Lund, Sweden 
Phone: + 46 - 46 - 222 48 16 
E-mail: torbjorn@thep.lu.se 
19936 + 11541 lines 



littp: / / thep.lu.se / tf2 / staff/ torbjorn / 



4.5.1 Introduction 

The Jetset program has been used frequently for QCD physics studies at LEP 1. For appli- 
cations at LEP 2, Jetset should be complemented with the Pythia program. While Jetset 
only gives access to one hard process, e"'"e~ 7*/Z'^ qq, Pythia contains a wealth of 
different processes. The two programs are fully integrated, in that a call to Pythia will not 
only generate a hard process but also automatically call Jetset routines to perform (timelike) 
parton showers and fragmentation. Output is in the normal LUJETS commonblock (with easy 
translation to the HEPEVT standard) and can be studied with the Jetset analysis routines. 
The emphasis of the Pythia/ Jetset package is to provide a realistic description of varying 
hadronic final states, but also non-hadronic processes may be generated. 

In addition to the briefer published description of the programs, there is a complete manual 
and physics description of over 300 pages |TB[. The programs, the manual, update notes and 
sample main programs can be picked up from the web address given above; additionally the 
CERN program library provides the programs and hardcopies of the manual. The description 
given here therefore only contains some highlights, with special emphasis on the aspects of 
relevance for LEP 2 applications. 

For the description of a typical high-energy event, a generator should contain a simulation 
of several physics aspects. If we try to follow the evolution of an event in some semblance of a 
time order, one may arrange these aspects as follows: 

1. Initially the e"^ and e~ are coming in towards each other. An electron contains virtual 
fluctuations into photons, quarks, gluons, and so on. It is useful to employ the same 
parton-distribution and parton-shower language as for hadrons. Thus also electrons and 
photons are included in the parton concept. An initial-state parton shower develops by 
branchings such as e — > e7, 7 qq and q qg. 
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2. One parton from each of the e"*"- and e~ -initiated showers enters the hard process, where 
then a number of outgoing partons/particles are produced. It is the nature of this process 
that determines the main characteristics of the event. (Also some soft processes are 
included in the program; since much of the same framework can be used we do not here 
belabour the differences.) 

3. If the hard process produces massive electroweak particles, such as the Z°, the or a 
Higgs, the decay into lighter objects must be considered. 

4. The outgoing partons may branch, to build up final-state showers. 

5. Further semihard interactions may occur between the other partons in the case of two 
incoming resolved photons. 

6. When a shower initiator is taken out of a beam particle, a beam remnant is left behind. 

7. The QCD confinement mechanism ensures that the outgoing quarks and gluons are not 
observable, but instead fragment to colour-neutral hadrons. 

8. Many of the produced hadrons are unstable and decay further. 

The time-order above does not have to coincide with the generation sequence. Typically the 
hard process is selected first. 



4.5.2 Hard processes 

Close to a hundred subprocess cross sections have been encoded in Pythia. Lepton, hadron 
and photon beams are allowed; thus the program can be used for pp/pp physics at the Tevatron 
or LHC or for ep physics at HERA. Here we concentrate on processes of relevance for LEP 2. 
Some of the more interesting ones are listed in table ^ and discussed below. Further comments 
may be found in other sections of this report. 

It it important to note that Pythia is not intended to be a precision program for electroweak 
physics. The philosophy is to provide sensible first approximations to a wide selection of hard 
processes, as a starting point for a detailed simulation of the subsequent QCD steps, i.e. parton 
showers, fragmentation and decay. It is therefore orthogonal in philosophy to many dedicated 
electroweak generators, that attempt to provide the hard-scattering cross section with very 
high precision but do not go beyond a parton-level description. 

Subprocess 1 is the familiar 7*/Z° process that dominates LEP 1 physics. The full in- 
terference structure between the 7 and Z° propagators is included. It supersedes the LUEEVT 
generator of Jetset. The main differences are: 

• LUEEVT uses a matrix-element approach to generate at most one initial-state photon, while 
Pythia allows for multiple photon emission in a parton-shower approach; 

• LUEEVT allows only hadronic final states, while Pythia also includes leptonic ones; 

• LUEEVT contains a simple Breit-Wigner with the width Pz as input, while Pythia contains 
an s-dependent Breit-Wigner that is dynamically calculated from electroweak parameters; 
and 

• the option to simulate first- or second-order MEs currently only exists with LUEEVT. 
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Table 8: Main LEP 2 physics processes available in Pythia. 



ISUB 



Process 



1 

18 
19 
22 
25 
35 
36 
69 
70 



10 
141 



Gauge boson production 



e^e 

e+e" 

e"'"e" 

e+e" 

e"'"e" 

67 - 

e7 - 

77 - 
7W 



^ 7VZ° 

77 

- 7(772°) 

^ (772°) (772°) 

^ W+W" 

e(772°) 



Higgs production 



24 


e^ 


^e 


103 


r 


T 


110 


e^ 




123 


e^ 


^e 


124 


e^ 
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Subprocess 19 contains a photon in addition to the 7*/Z . This means double counting, since 
already process 1 can contain initial-state-radiation photons, so results from the two processes 
should not be added. The usage of process 19 should be restricted to events that contain a 
high-px photon, where generation then is more efficient (and accurate) than what is offered by 
process 1. 

Subprocess 25 describes W pair production, including subsequent decay into four fermions 
with full angular correlations. The formalism includes s-dependent widths in the Breit-Wigners 
and options to pick the set of independent clcctroweak parameters. However, it is restricted to 
the basic graphs of W pair production ("CC03"). 

Subprocess 22 describes 7*/Z° pair production in a similar approximation ("NC02"). Note 
that interference terms between process 22 and 25 are not found anywhere. This is in accor- 
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dance with the basic philosophy of a reasonable but not exhaustive description of electroweak 
processes. 

Subprocesses 35 and 36 describe the production of a single 7*/Z° or W in the approximation 
of an effective photon flux. A process such as e"'"7 — > I^eW^ thus is convoluted with the parton- 
shower approximation of the e~ — > e~7 branching to give an effective process e^e~ — >■ Uee'W^. 
Process 35 has a singularity when the scattered electron has vanishing p±_ (in principle this 
is regularized by the electron mass, but in practice the rrie has been neglected). Therefore it 
is necessary to run with some minimum p± cut-off; for numerical reasons at least 0.01 GeV. 
An alternative description can be obtained by using an electron-inside-photon-inside-electron 
parton distribution (MSTP(12)=l) in process 1. For process 36 the decay of the W is assumed 
isotropic since the appropriate matrix elements have not been coded. Also subprocesses 69 and 
70 assume isotropic W/Z decay. Furthermore, process 70 does not include contributions from 
7* but only from Z°. The process implementations in this paragraph thus are less sophisticated 
than the single 7*/Z°, W pair and 7*/Z° pair processes above. 

Pythia is equipped with an extensive selection of production processes for the standard 
model Higgs, here denoted "h°". (It is called "H°" in the program, which confuses matters 
when two Higgs doublets are introduced, but for this report we stay with the conventional 
terminology.) Not all available processes have been listed, but only those of some interest. The 
most important by far (at LEP 2) is process 24, e+e" Z° ^ Z°h°. Both Z^'s in the graph 
have been included with a Breit-Wigner shape, so there is no formal restriction that either of 
them need be on or close to the mass shell. For a Higgs with + > E^m process 124 takes 
over, e"'"e~ — > ueVoW^W' — > uqUqIi^, but at a much smaller rate. Processes 110, 123 (Z^Z^ 
fusion) and 103 are even further suppressed. 

All major h° decay modes are included: h° ^ qq, h° h° ^ W+W~, h° Z°Z°, 

h° — > gg, h° 77 and h" — »• 7Z°. The branching ratios are automatically recalculated based on 
the Higgs mass. One point that should be noted is that the parton-shower algorithm matches 
to the same three-jet matrix element that is used for 7*/Z° decays. This gives a somewhat 
incorrect rate for three-jet production in Higgs decay. 

In the minimal supersymmetric extension to the standard model the number of production 
processes is further increased. The full set of Higgs particles is included in the program: h°, H°, 
A° and H^. Masses and couplings can be set by the user; this is a somewhat lengthy process, 
however, since currently the one-loop mass relations are not built into the program. The H° 
have the same production processes as the h°; the list in table ^ only shows the more interesting. 
Higgs pair production, H^H~, h°A° and H°A°, proceeds only through s-channel graphs and has 
been included as part of 7*/Z°/Z'° decays, process 141. The Z'" part can easily be switched off 
(MSTP(44)=4), so that processes 1 and 141 become identical except for the larger selection of 
decay modes in the latter. (Technically, this way the program can distinguish the Z' decaying 
to Higgses from a Z produced in Higgs decay, and accommodate different decay modes for the 
two.) 

Subprocess 141 is also useful for the study of virtual corrections caused by the existence of 
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a Z"' somewhere above the LEP 2 energy range. Vector and axial couphngs may be set freely 
to simulate various scenarios, and interference with 7*/Z° is automatically included. 

77 physics is a large area, in that a wealth of different subprocesses is involved. A photon 
may act as a pointlike particle or as a resolved, hadronlike state. A simple subdivision of 
processes is therefore into direct (58, 85), once-resolved (33, 54, 84) and twice-resolved (the 
rest, that is all processes allowed e.g. in pp collisions). The resolved part of the photon 
may be further subdivided into a VMD (vector meson dominance) and anomalous part. In 
total therefore six classes of events can be separated [|149|| . An automatic mix to provide a 
"minimum bias" sample of events is obtainable as an option (MSTP(14)=10). Processes 81-85 
include masses in the matrix elements, and thus are convenient to study e.g. heavy-flavour 
production. It should be noted that several aspects remain to be solved, for instance that of 
(slightly) off-shell incoming photons. Furthermore, on a technical note, Pythia is originally 
designed for fixed energies of the incoming particles, and so the process of having 77 "hadronic" 
collisions at varying energies is not yet fully automated. 

Finally, note that process 10 can be used both as a Bhabha and a deep-inelastic-scattering 
generator. In neither respect is it competitive with dedicated programs, but it may be useful 
for first estimates. 



4.5.3 Hard process generation 

The cross section for a process ij — > A; is given by 

aij^k = Jdxi J dx2ff{xi,Q'^)f] (x2, Q^) o-ij-^fc(s) . (13) 

Here a is the cross section for the hard partonic process, as codified in the matrix elements 
for each specific process. For processes with several particles in the final state it would be 
replaced by an integral over the allowed final-state phase space. The /j(x, Q^) are the parton 
distribution functions, which describe the probability to find a parton i inside an e^ beam 
particle, with parton i carrying a fraction x of the total e^ momentum, when the e^ is probed 
at some squared momentum scale that characterizes the hard process. The hard scattering 
therefore only involves a squared invariant mass s = X1X2S = XiX2E'^^, where E^m is the cm. 
energy of the event. 

The electron-inside-electron parton distributions are based on a next-to-leading order expo- 
nentiated description, see [Q. The approximate behaviour is 

f:ix,Q') - f (1 -x)i-^; /5 = ^ fln^ - iV (14) 

The form is divergent but integrable for x — 1, i.e. the electron likes to keep most of the 
energy. To handle the numerical precision problems for x very close to unity, the parton 
distribution is set, by hand, to zero for x > 0.999999, and is rescaled upwards in the range 
0.9999 < X < 0.999999, in such a way that the total area under the distribution is preserved. 
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In the 76 or 77 processes, an equivalent flow of photons is assumed, based on first-order 
formulae. There is some ambiguity in the choice of range over which emissions should be 
included. In the probably most appropriate alternative (MSTP(13)=2) the form is 



f;{^. Q') = ^ ^^^^ ""^^ In ( , ] . (15) 



e ' 



Here Q^ax (PARP(13)) is a user-defined cut for the range of scattered electron kinematics that 
is counted as photoproduction. Note that we now deal with two different scales, one related 
to the hard subprocess itself, which appears as the argument of the parton distribution, and 
the other related to the scattering of the electron, which is reflected in Q^^x- the default 
alternative (MSTP(13)=l) only one scale is assumed, i.e. Qmax(l ~ x)/x'^ is replaced by Q"^ 
above. 

Resolved photoproduction also involves the distributions of quarks and gluons inside the 



photon inside the electron. By default the SaS ID set [|15(]|] is used for the parton distributions 
of the photon, but several alternatives are available. 



4.5.4 Parton showers 

In every process that contains coloured and/or charged objects in the initial or final state, gluon 
and/or photon radiation may give large corrections to the overall topology of events. The philos- 
ophy of Pythia is to stay with the lowest-order cross sections (modulo trivial loop corrections 
such as the running of coupling constants) and then generate higher-order corrections in the 
parton-shower approach. This is less exact than the explicit calculation of higher-order matrix 
elements, but has the advantage that it can be applied also to processes where higher orders 
have not yet been calculated; additionally it includes multiple emissions. 

Showers may be subdivided into initial- and final-state ones, depending on whether they 
precede or follow the hard scattering. Of course, the subdivision often contains an element of 
arbitrariness, since interference terms may exist. In both initial- and final-state showers, the 
structure is given in terms of branchings a be, specifically e e'j, q qg, q qy, g 
gg, and g — i> qq. The kernel Pa-*bc{z) of a branching gives the probability distribution of the 
energy sharing, with daughter b taking a fraction z and daughter c the remaining 1 — z of the 
a energy. Once formed, the daughters b and c may branch in their turn, and so on. 

Each parton is characterized by some virtuality scale Q^, which gives an approximate sense 
of time ordering to the cascade. In the initial-state shower, spacelike values are gradually 
increasing as the hard scattering is approached, while timelike values are decreasing in the 
final-state showers. Shower evolution is cut off at some lower scale Qq, typically around 1 GeV 
for QCD branchings and around rrie for initial-state QED ones. From above, a maximum scale 
Qmax is introduced, where the showers are matched to the hard interaction itself. Unfortunately 
the selection of Qmax for a given hard scattering is not unique, but gives rise to some slop. 
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Despite a number of common traits, the initial- and final-state radiation machineries are in 
fact quite different. The Jetset final-state algorithm has been used extensively for Z° hadronic 
decays at LEP 1, and is not significantly altered since the LEP 1 writeup M. 



Initial-state radiation is handled within the backwards evolution scheme 111 5 111. In this 



approach, the choice of the hard scattering is based on the use of evolved parton distributions, 
which means that the inclusive effects of initial-state radiation are already included. What 
remains is therefore to construct the exclusive showers. This is done starting from the two 
incoming partons at the hard interaction, tracing the showers "backwards in time", back to 
the two shower initiators. In other words, given a parton b, one tries to find the parton a 
that branched into b. The evolution in the Monte Carlo is therefore in terms of a sequence 
of decreasing space-like virtualities and increasing momentum fractions x. Branchings on 
the two sides are interleaved in a common sequence of decreasing values. The definition of 
the X and z variables for off-mass- shell partons is not unique; in Pythia the z = Xb/xa of a. 
branching tells how much the scattering subsystem invariant mass-squared is reduced by the 
branching. If originally parton b was assumed to have vanishing p±, the reconstruction of the 
branching a ^ be introduces a p± for b, which is compensated by c. 



4.5.5 Beam remnants and multiple interactions 

The initial-state radiation algorithm reconstructs one shower initiator in each beam. Together 
the two initiators delineate an interaction subsystem, which contains all the partons that partic- 
ipate in the initial-state showers, the hard interaction, and the final-state showers. Left behind 
are two beam remnants. In some cases a remnant is a single object, as when a 7 is taken out 
of an e beam, leaving behind an e. When taking an e out of an e, a soft 7 is left behind, which 
is then more related to the cutoff of Q^) at x = 0.999999 than to the ordinary beam- 
remnant concept, but is handled with the same machinery. In other cases a remnant consists 
of two objects, as when a q is taken out of an e, leaving behind e -|- q. The latter example has 
a coloured remnant, meaning that the fragmentation of the hard-process partons is connected 
with that of the beam remnants. 

A resolved photon contains many partons. In a twice-resolved 77 event there is thus the 
possibility of multiple interactions, i.e. of multiple semi-hard parton-parton processes in the 
same event. A model for this phenomenon is included in Pythia [p.49|| , and will be further 



developed to better represent differences between the VMD and anomalous states. 



4.5.6 Fragmentation and decay 

The Lund string fragmentation description |^ and the decay routines in Jetset have not 
changed significantly since the LEP 1 writeup 0, and so are not described here. The string 
fragmentation approach has been generally successful in comparisons with LEP data, although 
some shortcomings have shown up. See section || for further details. 
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The issue of Bose-Einstein effects has received increased attention in recent years, e.g. in 



connection with possible consequences for the W mass determinations [|152|| . The existing al- 
gorithm 1^ works well in many respects, but is by no means to be considered as a definite 
solution to the problem. A somewhat different approach has been implemented to allow some 
cross-checks, and further alternatives may appear in the future. In the current standard algo- 
rithm, identical particles are pulled closer together in such a way as to enhance the two-particle 
correlation at small relative momentum separation. This makes jets slightly narrower, so that 
fragmentation parameters have to be retuned for reasonable agreement with data. In the al- 
ternative, the shift of identical particles is somewhat reduced, while non-identical particles are 
pushed apart a bit, so that the average properties of jets remain unchanged. This alterna- 
tive does not yet come with Jetset, but is available as a plug-in replacement for the LUBOEI 
routine, at |http:/ /thep.lu.se/tf2/stafl:"/torbjorn/test/mainl0.1| . 



Also colour rearrangement has been extensively discussed in recent years. Code that allows 
this has not yet been integrated in the standard Pythia/ Jetset libraries, but is obtainable 



separately, see section 4.7 



4.5.7 Final comments 

Pythia/ Jetset are likely to be among the major event generators at LEP 2: access to a 
broad selection of hard scattering subprocesses is combined with a well-tested description of 
parton showers and fragmentation. Limitations exist, however. Pythia is not a program for 
precision extraction of electroweak parameters; for instance, no (non-trivial) loop corrections are 
included in the matrix elements. One may well imagine hybrid arrangements, where dedicated 
generators are used to provide an improved description of some especially interesting hard 
scattering processes, such as four-fermion final states, while the rest of the Pythia/ Jetset 
machinery is used to turn a simple parton configuration into a complex hadronic final state. An 
example of such an interface is discussed in section |5.3| . Furthermore, the Ariadne program 
for colour dipole radiation offers an alternative to the parton-shower description of Jetset, 
and can be used for all the hard processes in Pythia. 

While there are no major additions planned for Pythia/ Jetset, the intention is to continue 
a steady development and support activity. The 77 sector maybe is the area where most further 
studies are required to complete the picture, but also other aspects deserve attention. 



4.6 UCLA ansatz 



Basic Facts 

Program names: 

Versions: 

Author: 



UCLA m |T53| 



UCLA 7.41 of 1 October 1995 
Sebong Chun and CD. Buchanan 
Department of Physics 
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UCLA 

405 Hilgard Ave, Los Angeles, CA 90024 
USA 

Phone: (310) 815-1992, 7466 
E-mail: chun@physics.ucla.edu, 

buchanan@physics.ucla.edu 
Program size: 1922 lines 

Program location: [http: / /www. physics. ucla.edu/ ~chuns 



The goal of the UCLA hadronization modeling is to study and develop the underlying 
principles of e"'"e~ annihilation into hadrons, constructing a simple phenomenology which can 
be used both as a "target" for non-perturbative QCD calculations and also to accurately predict 
data. 

The UCLA7.41 program, a spin-off of the Lund relativistic string Monte Carlo program 
Jetset, is the manifestation of this modelling to be used in comparing predictions with e^e~ 
data. As Jetset has upgraded to new versions, the UCLA program has likewise been adapted 
with a parallel nomenclature. 



The modern UCLA modeling presumes that, by making a few assumptions which can 
be rationalized within a QCD context (for example, a strong coupling expansion in lattice 
QCD), one can construct a Weight Function for any specified e+e~ — > hadrons event. That is, 
given the center-of-mass energy of the e+e~ system and the flavor and momenta of the primary 
hadrons produced, the UCLA modeling attaches a weight to the entire event, to be used in 
comparison with other possible events at that Ecm- 

The general structure of the Weight Function (in addition to kinematics of energy/momen- 
tum conservation and phase space with limited transverse momentum) depends on (a) an area 
law in space-time, (b) possible suppression factors at the vertices where a virtual qq pair is 
created from the colorfield, (c) "knitting factors" to knit a quark and antiquark together into 
the spatial wave function of a meson (or quark and diquark into a baryon), and (d) Clebsch- 
Gordon coefficients to knit the quark and antiquark (diquark) together into the flavor and spin 
state of the meson (baryon). 

a) The area law is exp(— 6'y4) where b' is a constant and A is the area enclosed by the 
quark and antiquark trajectories in a space-time plot of the event. Almost any strong-coupling 
interaction will, in fact, give this sort of dependence. 

b) The UCLA modeling assumes that there is no signiflcant vertex suppression for qq pairs 
if the quark mass is less than the hadronic scale of ~ 1 GeV; that is, uu, dd and ss all have 
probability ^ 1.0 of virtual creation from the colorfleld. 

c) The UCLA modeling assumes that all knitting factors are comparable, whether the 
hadron to be constructed is a spin or 1 meson or a spin 1/2 or 3/2 baryon. (Probability 
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normalization on the fragmentation function derived below yields a value of the knitting factor 
of ~ (40 Mev)-2.) 

d) The Clebsch-Gordon coefficients are simply the relevant flavor/spin coupling of a quark 
and antiquark (diquark) into a meson (baryon). Note that (c) and (d) taken together describe 
the coupling of a quark and antiquark or diquark into the complete state function of a hadron. 

Although, in principle, knowing the Weight Function for the final state is enough to select an 
event, it is practically impossible to implement in this form. In order to implement this simple 
event Weight Function approach into a working Monte Carlo program, it is necessary to derive 
a fragmentation function for an "outside-in iterative one-particle-at-a-time" implementation 
such as Jetset uses. 

By somewhat lengthy but straightforward algebra, this can be accomplished. The result 
so derived turns out to be the Lund Symmetric Fragmentation Function (LSFF) [|153|| , with 
normalizing parameters of the vertex suppression, the spatial knitting factor, and the Clebsch- 
Gordon coefficient (see |^^). That is, the UCLA modeling simply amounts to using the LSFF as 
a hadronic production density weighted by Clebsch-Gordon coefficients, where the suppression 
of heavy mass particle production arises entirely from the exp(— 6m^/z) factor in the LSFF. 
(Note: the general structure of the Weight Function and the subsequent derivation of the 
fragmentation function can also be used to describe the Lund Jetset treatment. The difference 
is that Jetset presumes an ss vertex suppression of about 0.3 and a knitting factor for vector 
mesons of about 30% of that for pseudoscalar mesons, does not in general use Clebsch-Gordon 
coefficients, and adopts a normalization scheme that does not incorporate the exp{—bm'^/z) 
factor.) 

The UCLA7.41 program uses the parton shower and decay table parts of Jetset, but 
replaces the flavor and momentum selection part with the UCLA modeling ansatz described 
above. Default values for the parton shower are A = 0.2 GeV and Qo = 1.0 GeV. Meson 
production is controlled by the two natural parameters of the LSFF with default values of a = 
2.1 and b = 1.1 GeV~^. Local transverse momentum compensation is approximated by a factor 
of exp{— ^^^bp^/ z) , where n is a parameter of default value 2.0. For baryon production, with 
"popcorn" mesons produced between baryon and antibaryon, an additional popcorn suppression 
factor of exp(— ?7mpop) is introduced with the default value of r/ = 10 GeV~^. For more details, 
please refer to refs. 



This structure and values gives a rather good description of multiplicities, inclusive distri- 
butions, and correlations for hadron production from E^m of 10 to 91 GeV, with the possible 
exception of the spin 3/2 baryons at 91 GeV. The description of heavy flavor (c, b) production 
distributions also seems reasonably good, with no additional parameterization or parameters. 

For a detailed instruction on how to set up parameters and use the program, please re- 
fer to the manual at WWW location |http: / / www. physics. ucla.edu/^chuns| . A short set of 



instructions is available in the header to the actual program. 
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4.7 Colour reconnection codes 



One of the QCD questions that has attracted attention in recent years is that of colour recon- 
nection (or colour rearrangement) ||127| , |154| , p.55| , |156|| . This issue has implications for W mass 



studies, but is also of interest for our general understanding of QCD. 

The concept may be illustrated by the process e"'"e~ W^W~ — > qiq2 q3q4. To first ap- 
proximation, the hadronic final state can be viewed as coming from the incoherent superposition 
of two sources of particle production: the q^cJa and the q3q4 ones. If colours are reconnected, 
the sources would instead be q^q^^ and q3q2- The picture is complicated by the possibility of 
gluon emission. Gluons with an energy above the W width can be viewed as independently 
emitted from the respective W source, to a good first approximation: propagator effects ensure 
that interference terms are suppressed. No similar suppression exist for soft gluons or in the 
nonperturbative regime. Therefore standard calculational techniques are of limited interest, 
and the phenomenon mainly has to be studied within the context of specific models. By now, 
several independent codes exist, some part of existing QCD generators, others available as 
add-ons. Below we list the known ones and give some specific details. 



4.7.1 A PYTHIA-based implementation 

(code by T. Sjostrand) 



The code used for the studies in ||155|| has not (yet) been incorporated in the Pythia/ Jetset 
programs. A sample main program and the colour rearrangement subroutines can be obtained 
at web address [littp: / / thep.lu.se / tf2 / staff'/ tor bjorn/ test / mainOl.l] . Several different options are 



available, among others: 

• scenario I, where strings are considered as extended colour ffux tubes and the reconnection 
probability is proportional (up to saturation corrections) to the space-time overlap of the 
W"*" and W~ strings; 

• scenario II, where strings are considered as thin vortex lines and reconnection may occur 
when strings cross; 

• scenario II', a variant of scenario II where only those reconnections are allowed that reduce 
the total string length; and 

• the instantaneous scenario, where reconnections are allowed before the parton-shower 
evolution ||154|] ; unphysical but handy for comparisons. 



At most one reconnection is performed per event, in scenario II the one that occurs first in time, 
in scenario I selected according to relative probabilities given by the overlaps. Reconnections 
within a W system are not considered. 



4.7.2 Another PYTHIA-based implementation 

(code by S. Todorova) 



The code follows the physical approach of ||155|| . The reconnection phenomenon is simulated 
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with the help of the string model, where strings are considered to be either flux tubes or vortex 
lines with arbitrary diameter of the core. In order to get a more realistic estimation of the 
effect of colour reconnection, the following features (not found in the preceding code) were 
incorporated in the simulation: 

• space-time evolution of parton shower; 

• multiple reconnections; and 

• self- interaction of a string (production of glueballs). 

It should be noticed that the space-time evolution of a shower together with the self-interaction 
of a string allow the study of string reconnection effects in a single parton shower (decay of a 
single TP). 

The search for candidates for reconnection is processed in parallel with the shower devel- 
opment (reconnection can take place before the emission of the last partons). Overlaps of the 
colour fields of flux tubes are calculated numerically (using multichannel MC integration with 
importance sampling). The method is slow but this is the price to be paid for (relative) accu- 
racy and individual treatment of each event. The minimal distance between "vortex lines" is 
found by a minimization procedure based on parabolic fit. 

The code is available in the directory crnvax:[nova.colour_ reconnection]. 



4.7.3 An ARIADNE-based implementation 

(code by J. Hakkinen) 



The aim of the simulation program presented in [|156|| is not so much to study the effect 
of recoupling on the average events, but to study if rare recoupled events can be identified. 
Perturbative QCD favours states which correspond to "short strings" , i.e. parton states which 
produce few hadrons. This string "length" can here be specified by the A measure, defined 

m 



m 



which correspond to an effective rapidity range. If recoupling occurs it is conceivable 
that it is favoured when the recoupling produces a state with lower A measure, and such states 
may also be more easy to identify. For this reason the program produces recoupling such 
that the A measure for the reconnected final state is minimized. Gluons with E ^ F-yv are 

158 , 155 1 . This emission is simulated 



For gluons with cm. 



emitted independently within the original qq systems 
using the Dipole Cascade Model [|122|] implemented in Ariadne |13 
energy below Fw ~ 2 GeV there may be unknown interference effects due to emission from 
the two W systems. These low-energy gluons give very little effect on the hadronic final state, 
however, if the hadronization phase is described by the Lund string model |]28| implemented in 
Jetset [15|. They are therefore disregarded in the parton states, which implies that a small 
fraction of the energy {k. 4%) will be lost in the simulations. During minimization of A all 
possible final state configurations, obtained by cutting the original gluon chains in one place 
only before reconnecting to two new systems, are compared with each other. In this program 
reconnection between the two strings occur once in every event while reconnections within the 
strings are not considered. Thus, the program is not expected to reproduce average events, but 
possibly a small admixture of recoupled events. 
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The C code used in ||156|| is available at [http : / / 1 hep .lu.se/tf2/ hep / hep . litmll or through 
anonymous ftp at thep.lu.se:/pub/LundPrograms/Misc/wwpair.tar.Z. The code contains two 
more models; the instantaneous scenario of | |154|| , and random reconnection of the strings. 
These models are only used for comparison with the "main" model. 



4.7.4 Another ARIADNE-based implementation 

(code by L. Lonnblad) 



The model in ||127|| for colour reconnections, implemented in the Ariadne program 



is similar to the one in ||156|| in that it reconnects colour dipoles within the framework of the 
Dipole Cascade Model (DCM) ||122|| with a probability l/N^, only if the total string length 



becomes reduced. The main differences are that reconnections within each W system (and also 
among the partons from a Z decay) is allowed, that several such reconnections are allowed in 
each event, and that reconnections are allowed during the perturbative cascade. 

To achieve this, colour indices are assigned to each dipole, and after each emission, dipoles 
with identical indices are allowed to reconnect. The indices are chosen randomly, but restrictions 
are made to ensure physical colour flows, e.g. two gluons created by a gluon splitting should 
not be allowed to form a colour singlet. In the DCM, however, a gluon is radiated coherently 
by the dipole between two partons, and a procedure has to be introduced, where the emitted 
gluon is said to have been radiated off one of the two emitting partons with some probability 
depending on which is closer in phase space. 

In the case of e"'"e~ W^W^ reconnections are initially only allowed within each W 
system separately. After all gluons with > Fw have been emitted, reconnections between 
the W systems is switched on and gluon emission with Eg < Fw is performed in the possibly 
reconnected systems before hadronization. 



4.7.5 A HERWIG-based implementation 

(code by B.R. Webber) 

A model for colour reconnection has been implemented in a package of subroutines that 
can be used with HERWIG (version 5.8). The new integer parameter IREC0=0,1,2 determines 
the reconnection option used. IRECO=0 means no reconnection and IREC0=2 gives "immediate" 
reconnection of the quark-ant iquark pairs in hadronic WW events, before parton shower gener- 
ation, with probability PRECO (default value = 1/9). In HERWIG this changes the evolution of 
the showers, as well as the colour connections, because the initial opening angles are different. 

The most serious option is IREC0=1, which invokes a model based on the assumption that 
reconnection occurs locally in space-time. First, a space-time structure is computed for each 
parton shower in the event. This is done using a package written by Mike Seymour to store 
the internal lines of showers, which is turned on by setting INTLIN=.TRUE. . The algorithm is 
semi-classical, but qualitative features and orders of magnitude should be correct. In the case 
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of W hadronic decays, each W decay point is generated with the appropriate exponential decay 
distribution. Then the locations of all the vertices in the showers are computed by assigning a 
space-time separation Axj — qi/ {qf — m]) to vertices joined by an internal line i of 4-momentum 
qi and on-shell mass rrii. 

The HERWIG cluster hadronization model, normally called immediately after showering has 
terminated, involves splitting each final-state gluon into a quark-antiquark pair. For each quark 
i there is a colour partner antiquark j, with which the quark would normally be paired to form 
a colour singlet cluster (ij). The IREC0=1 option introduces a reconnection phase before cluster 
formation. In this phase the program looks for another colour-connected quark-antiquark pair 
{kl) such that (il) and {kj) could be colour singlets and 

\Axii\^ + \Axkj\^ < \Axij\^ + \Axm\^ , (16) 

where Axij is the (ij) cluster size, defined as the separation of the production vertices of i and 
j (note that this can be zero, e.g. if i and j come from a W decay that did not radiate any 
gluons). If such a pair exists, the reconnection {ij){kl) {il){kj) would reduce the cluster 
sizes, and so it is performed with probability PRECO. Note that reconnection can happen inside 
a single shower and not just between different showers. Thus some retuning of parameters to 
fit data on e+e~ ^ Z° — > hadrons will be necessary when using IREC0=1. 

The code can be obtained by anonymous ftp from 
hep.phy.cam.ac.uk = 131.111.66.27 
The following files should be copied from directory disk$alphal: [public .herwig] : 

• hwwmas58.for - sample main program and analysis routines; 

• hwreco58 .for - modified HERWIG routines HWBFIN, HWBJCO, HWCFOR which replace those 
in HERWIG version 5.8, plus new routines HWGCLU, HWGCMO, HWUPIP, HWVHEP. 

There is a new common block containing relevant parameters and counters: 

COMMON/HWRECO/PRECO , EXAG , IRECO , MEVTS , MCLUS , MRECO , MSWCH , INTLIN 

PRECO is the reconnection probability (default 1/9); EXAG is an 'exaggeration factor' for the 
W lifetime, to study effects of the WW separation (default 1.0); IRECO is the reconnection 
option (see above, default 1); MEVTS etc. are integer counters for number of events, clusters, 
reconnections, and WW reconnections; INTLIN is set to .TRUE, when IREC0=1 (see above). 

The code is still under development; please notify webberQhep . phy . cam .ac.uk of any prob- 
lems, bugs and/or peculiarities. 
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4.8 Monte Carlo Implementations of Exact Next-to-Leading Order 
Calculations 

Basic Facts 

Program name: EVENT EERAD EVENT2 

Authors: Zoltan Kunszt Walter Giele Stefano Catani 

Paolo Nason Nigel Glover Mike Seymour 

email: nason@suryall.cern.cli E.W.N.Glover@dur.ac.uk seymour@suryall.cern.cli 



There are now three publicly- available programs for calculating next-to-leading order cor- 
rections to arbitrary infrarcd-safc two- and three-jet quantities in e^c^ annihilation. Although 
these use Monte Carlo integration techniques, they should be contrasted with Monte Carlo 
Event Generators in several ways. Firstly, they calculate the exact result in perturbation the- 
ory for the C(q;s) corrections to a given quantity — no more nor less. Secondly, the phase-space 
configurations generated do not have positive-definite weights, so a probabilistic interpretation 
is not possible. Finally, for both these reasons, the programs only ever consider the partonic 
final state, and no treatment of hadronization is attempted. 

There arc many advantages of implementing higher-order QCD calculations as matrix- 
element Monte Carlo programs. For all but the simplest observables, the required phase-space 
integrals are not analytically tractable, and some form of numerical integration becomes manda- 
tory. Since each phase-space point sampled by the program has a direct correspondence to a 
set of final-state momenta, any infrared-safe jet or event-shape definition may be used, and can 
be implemented exactly as in an experimental analysis. Many event properties can be analyzed 
simultaneously, simply by adding code to the analysis routine of the program to histogram the 
quantity of interest. 

However as is well-known, the real and virtual corrections are separately divergent but 
with finite sum, so naive numerical integration of each matrix element would fail. Thus a 
regularization scheme must be used to render the integrals finite. It is principally in the 
definitions of regularization scheme that the three programs differ, although there are other 
important differences. 

The diff'erence between the regularization schemes can be illustrated using a simple one- 
dimensional example. In dimensional regularization using (4 — 2e) dimensions, the integrals we 
encounter are typically of the form 

where the first part represents the real cross-section, the second the virtual, and x would 

typically be a gluon energy or parton-parton invariant mass. The function 0{x) represents a 
final-state observable that is infrared safe, i.e. with the requirement that 0{x) tends smoothly 
to 0(0) as X tends to 0. If the integral were analytically tractable, it would yield an e-pole 
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that canceled the virtual term, leaving a finite result. However, in general it is not, and we 
must manipulate it into a form in which the physical limit e — can be taken before numerical 
integration, without making any assumption about 0{x). 



The phase-space slicing method does this by introducing an unphysical parameter xq, 

'Xq X t Jxq X 



(0)=/ ^p{x)+ O(x) + -O(0)^/ — 0(x) + log(xo)0(0) + 0(xoO'(0)). (1^ 

JO X Jxn X e Jxn X 



The result becomes exact in the limit Xq — 0, which practically means Xq <^ a^physicai, where 
^^physicai is the smallest physical scale in the problem. 

The subtraction method works by subtracting and adding a term derived by projecting each 
point in four-parton phase-space onto some point in three-parton phase-space, and calculating 
the observable at this phase-space point together with an approximate matrix element. This 
must be such that it matches all the divergent terms of the full matrix element. In our simple 
example this corresponds to 

(O) = /' 4^0(x) - /' 4^0(0) + /' 4^0(0) + -O(O) = /' - (0(x) - 0(0)) . (19) 
JO x^+^ JO x^+^ JO x^+^ e Jo X 

Note that this is exact and does not depend on any unphysical parameters. 

The matrix elements for 7* — qqg have been known to next-to-leading order for many 
years ||15^^] . These were later checked by other groups, and used for specific calculations of a 
variety of event shapes. For the 'QCD at LEP' report ||16CI|| , Kunszt and Nason wrote a general- 
purpose Monte Carlo program using the subtraction method that could calculate the next-to- 
leading correction to any event shape or jet definition, EVENT[|161||. This has been considered 



the standard calculation for many years, but has two significant shortcomings owing to the 
matrix elements used: they have been summed over permutations of the outgoing partons, which 
means that quarks and gluons cannot be distinguished in the final state; and they consider the 
decay of a virtual photon, so can only predict quantities averaged over orientations of hadronic 
events, losing all information on their lab-frame directions and lepton-hadron correlations. 
Furthermore they neglect specific axial-axial contributions that as a point of principle are 
essential for describing Z° decays, although in practice these are never numerically significant. 

More recently two groups have proposed general algorithms for calculating next-to-leading 
order corrections in arbitrary processes, and both have used three-jet production in e"'"e~ anni- 
hilation as a simple first proving ground for their methods. These have resulted in the EERAD 
program by Giele and Glover [^6^, which uses the slicing method, and the imaginatively-titled 
EVENT2 program by Catani and Seymour [[L63|] , which uses the subtraction method. Both of 
these use the full next-to-leading order matrix elements for e"'"e~ — >■ qqg, avoiding the short- 
comings of EVENT. Although EVENT2 uses the matrix elements of the Leiden group |ll64| by 



default, it has options to use the same matrix elements as EVENT or EERAD as a cross-check. 



Numerical results of the three algorithms are discussed in ||165|| and shown to be in good agree- 
ment. Since they are supposed to be exact calculations of the same quantity, rather than 
models, any differences between them should be treated as bugs. 
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5 Standardization 



5.1 Particle codes and /HEPEVT/ update 



The /HEPEVT/ standard ||166| | has been widely adopted by Monte Carlo authors for storing 



information on generated events. In practice the real variables are commonly declared to be 
DOUBLE PRECISION and often the size is expanded to NMXHEP=4000. We propose that these are 
now added to the standard. 

In /HEPEVT/ it was intended for particles to be identified using the PDG numbering scheme 
[|167|| . However the conventional numbers assigned have deficiencies, particularly concerning the 



neglect of particles expected according to the quark model but not yet identified in experiment, 
for example the hb. This proves troublesome for those program authors who include such states 
and has lead to ad hoc solutions. Further the higher, orbitally excited L = 2, 3, . . . and radially 
excited n = 2, 3, . . . mesons are labelled in a somewhat unsystematic way. In order to preserve 
the concept of uniqueness, allow for the missing quark model states, systematize the numbering 
and remain true to the spirit of the PDG scheme we suggest the following revised numbering. 

Table ^ lists the n = 1, L = 0, 1 mesons and indicates their numbering, for these states 
this is largely in accord with the PDG scheme and with the stdhep (Jetset) implementation 



168|| . In the pairs of / = 0, (u, d, s) mesons: {r],r]'), (a;,0), (hi(1170), hi(1380)), etc. the 
lighter state is labelled 22 and the heavier 33, refiecting the naive, dominant quark contents. 
Bound states involving top quarks are not expected, due to the quark's high mass, and therefore 
are not considered. The mixed Kg and K^^ states are still labelled 310 and 130 respectively. 
The table should be extended to include n = 1, L = 2,3, . . . states; this leads to up to four 
mesons of the same total spin. It is proposed to reserve the fifth digit to differentiate these 
states by continuing the sequence established for the L = 1 mesons. That is, for a given J > 
the numbers would be: {L, S) = {J — 1,1) : -k -k m, (J, 0) : 10 -k km, (J, 1) : 20 * km and 
( J + 1, 1) : 30 k km, where as usual m = 2 J + 1. The J = states represent an exceptional 
case, here we propose {L, S) = (0,0): kkl and (1,1) : 10 k kl, as done in table ^; this 
may be thought of as LO -k-kl. Radially excited mesons, n = 2,3..., are effectively copies of 
the the above states, it is proposed to introduce a sixth digit to differentiate them as follows: 
n = 1 : k-kk-kk, n = 2 : l-k-k-k-kk, n = 3 : 2 -k-k-k-k-k, etc. Thus for example the K*~^ (1680), a 
l^Di state would be numbered 30323 and the p°(1450), a 2^5*1 state 100113. The numbering 
of excited mesons suggested here differs significantly from the original PDG scheme. 

Table ^ lists the lowest lying J = 1/2, 3/2 baryons, including the anticipated charm and 
bottom states. Two J = 1/2 states exist for baryons containing three different fiavours of 
quarks. When the two lighter fiavours are in a symmetrical (J = 1) state the baryon is called 
a S, S' or Q' and a A, S, or Q if they are in an antisymmetric (J = 0) state. To distinguish 
the lighter, antisymmetric states the light quark numbers are reversed; for example^ has 

''Actually the PDG naming rules do not make it clear which state to put the prime on, we have provisionally 
chosen to place it on the heavier state. 
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Table 9: Proposed numbering scheme for the lowest lying mesons: for example the af has the 
number —20213. The names of the pseudovector particles are distinguished by their masses, if 
these are not currently estabhshed L and H are used to indicate hght and heavy. The pseudo- 
vector states Ki(1270) and Ki(1400) are beheved to be admixtures of the Ki^ l^Pi and Kia 
l^Pi, in such situations the lighter state is given the lower number. 

number 5312 and the number 5132. This extends the convention that heavier states are 
given larger numbers. Excited states are not yet incorporated into event generators and thus 
are not covered here. 

Increasingly supersymmetric particles are found in event generators, we therefore take his 
opportunity to put forward the following numbering scheme for them. A seventh digit is added 
being: either 1 [1 -k -k -k -k -k -k) for the partner of a boson or left-handed fermion; or 2 [1 -k -k -k -k -k -k) 
for the partner of a right-handed fermion. When left-right mixing occurs the ordering should 
be by mass. Examples include: 
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Table 10: The proposed numbering scheme for the baryons. In the first J = 1/2 column the 
order of the light quark numbers is reversed; for example A has number 3122 whilst E° is 3212. 



The possibility of numbering potential SUSY mesons and baryons in the same spirit is left open 
at present. 



5.2 Decay Tables 

The study of identified particle production and the physics Tinderlying the hadronization mech- 
anism continues to be an active area of research at LEP 1. In hadronic Monte Carlo event 
generators final state particles are produced in two stages. Primary hadrons come directly from 
the clusters/strings/ e^c. that model the non-perturbative parton to hadron transition. Subse- 
quently chains of secondary particles arise from the decays of the unstable primary hadrons. 
Note this separation is well defined in the context of MC programs but in reality for the short 
lived, strongly decaying resonances it may be only semantics. 

Thus a major common component of hadronic MCs are routines to do the decay of unstable 
particles. These are based on the use of tabulated branching ratios and basic matrix elements, 
though in the case of r's and b-hadrons speciahzed packages are also available. The construction 
of these decay tables is not a simple task and requires much per-in-spiration to fill in gaps in 



65 



present measurements ||169|| and deal with problematic cases. It would save much duplication 
of effort if one basic table could be used by all programs. This implies the ability to swap decay 
tables and thereby would allow some control over a (spurious) source of apparent variation in 
the rates of primary hadron production in the different hadronization models. A common, user 
friendly, interface would also enable easy maintenance and modification of the tables by users. 

To achieve such a goal requires a unique way of identifying the particles, and any associated 
matrix elements, together with a standard format for outputting and inputting the tables. The 
revised PDG codes above provide a unique and logical means of identifying the particles. To 
identify the matrix elements we propose developing a set of standard three-digit integer codes. 



following the convention of table 11. 



Code 


Matrix Element 





Isotropic decay 


1-99 


Standard codes to be agreed 


> 100 


Program specific options 



Table 11: Proposed convention for matrix element codes 



It is reasonable to restrict both the number of decay products to five, using zeros to complete 
an entry, and also numerical branching ratios to five decimal places. A more than five body 
decay can be stored, realistically, as a sequence of decays involving intermediate resonances. In 
studies involving very rare decays it is sensible to use a higher branching ratio and then apply a 
compensating normalization factor. It is then proposed to write out the following information. 

Number of decays listed 

Decaying particle, branching ratio, matrix element code, 1-5 decay products 
using the following FORMAT statements, 
100 FORMAT (IX, 14) 

200 F0RMAT(1X,I8,1X,F7.5,1X,I3,5(1X,I8)) 

An example is provided by the 7r° decays: 
2 

111 0.98800 22 22 ^ 77) 

111 0.01200 101 22 11 -11 (7r0^7e-e+) 

It must be recognised that b-hadrons represent a special case. In the absence of detailed 
knowledge about a significant fraction of their decays MC programs resort to models based 
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on partonic decays and fragmentation. Partonic decay modes may also be stored in the above 
format. Suppose in a bq hadron (here q may represent a diquark) the decay is b ^ cW~ — >■ 
cqiq25 this can be coded in one of two sequences either 'qi,q2,c' or 'c,q2,qi'. These two 
options can be exploited to refer to the two possible colour connections separately: (cq)(qiq2) 
and (cq2)(qiq) respectively, at the discretion of program authors^. 

It is now simply a matter of providing a .DAT file containing the decay table listed in the 
above format. To standardize the interface to the individual MC programs the following two 
subroutines are proposed: 

7^7^I0DK(IUNIT,IF0RMAT,I0PT) and 
T^MODK ( IDK , BR , ME , IPRDl , IPRD2 , IPRD3 , IPRD4 , IPRD5) 

were identifies the MC program. The first is used to read lUNIKO or write IUNIT>0 the 
decay table to the given unit number with IFORMAT specifying how the particles are identified. 
The standard is IF0RMAT=1, that is use the revised PDG codes; nonportable program specific 
options may include: =2 use the internal numbering or =3 use the internal character string 
names. Authors and users may prefer the later options as more transparent than the PDG 
numbers. If I0PT=1 then matrix element codes > 100 (program specific) are accepted, if 
IOPT=0 then such codes are treated as not recognised and set to zero, isotropic decay. The 
subroutine *T*rMODK is intended to allow individual lines of the table to be modified or added, 
before or during event generation; the arguments follow the standard format. Note that when a 
new mode is added or an existing branching ratio modified the sum of the remaining branching 
ratios should be rescaled to preserve unit sum. This means that when two modes of the same 
particle are altered the order of the calls is important for their resultant branching ratios. 

The provision of such an interface rests with the actual program authors who need to 
convert between the standard format and their own internal structures. These interfaces may 
be expected to be robust against unrecognized or blank particle names and provide basic checks 
of the allowed kinematics, electric charge conservation and unit sum of branching ratios. Such 
an interface has been established for HERWIG and successfully used to import the Jetset 
decay table. However if program users do modify the provided decay tables then they must 
accept responsibility for them making sense. 

5.3 Interfaces to elect roweak generators 

A number of dedicated four-fermion generators are being written for LEP 2 applications. The 
good ones will do the electroweak theory much better than standard general-purpose QCD 
generators. On the other hand, they do not contain any QCD physics aspects, i.e. neither 
perturbative parton showers nor nonperturbative hadronization. This makes the electroweak 

^Observe that if the V — A matrix clement was constructed as (po-P2)(pi-P3) these would both give 
(b.q2)(qi-c). 
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(EW) generators well suited for some applications, such as total cross sections and leptonic final 
states, but generally unsuited for the study of hadronic or mixed hadronic-leptonic final states. 
It is therefore logical to interface them with parton-shower and hadronization programs. To 
some extent, this is already happening. However, in writing these interfaces there are certain 
dangers involved. There may also be a lot of work involved. 

It would therefore be advantageous if the event generator authors involved could agree on 
a common approach: EW authors provide the four-fermion configuration in a standard format 
and QCD authors provide a standard interface that converts this to a set of final hadrons. Then 
only one interface needs to be written for each program, instead of one for each combination of 
EW and QCD programs. In this section we propose such a standard and report on progress in 
implementing it. 

5.3.1 The basic problem 

In the electroweak sector, fermions can be viewed as asymptotically free final state particles. 
This means that the production of a specific final state is fully calculable perturbatively. Many 
different intermediate states can contribute to the same final state, without any ambiguities 
being introduced by this. The total probability for a final state is given by squaring the sum 
of amplitudes 

|A|2 = |A+A2 + ... + ^nr (20) 

(suppressing the issue of helicity sums, etc. — these aspects are not important for the general 
discussion) . Interference effects therefore are included automatically. 

QCD is different. Quarks are not asymptotic states. The final state consists of colour singlet 
hadrons, not coloured partons. The transition from perturbative to non-perturbative physics is 

not understood from first principles, but is at present modelled. The model used describe well 
what happens to a simple quark- ant iquark pair, e.g. — > qq at LEP 1. At LEP 2, four-quark 
states qiq2q3q4 have to be mastered. If we want to make use of our hard-won phenomenological 
experience, it is therefore essential that the qiq2q3q4 system can be subdivided into two colour 
singlet subsystems, either q^qg -|- €[3^4 or q^q4 -|- q3q2. Each subsystem can then be described in 
the same way as a LEP 1 event. On the contrary, if we are not allowed to use such a subdivision 
into singlets, a completely new hadronization formalism would have to be invented (with brand 
new parameters to be tuned to the LEP 2 data). 

Unfortunately, there are complications. As a simple illustration, consider a system uddu. 
The production obtains contributions from several possible intermediate states. One is a W^W~ 
pair, with W"*" — > ud and W~ — > du. Another is a Z°Z° pair, with the first Z° — > uu and the 
second — > dd. These two alternative intermediate states correspond to different colour 
singlets, and therefore would differ with respect to the treatment of subsequent parton showers 
and hadronization. That is, the final state contains a "memory" of the intermediate state. 
Furthermore, \A\'^ contains an interference term between the two alternatives, where the colour 
flow is not well-deflned in perturbation theory. This is reflected in a relative colour factor 
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l/(iV^ — 1) for the interference term, meaning e.g. that both uu and ud are in relative colour 
singlet states. Kinematical factors are not likely to compensate for the colour suppression, 
so numerically the interference terms may not be large. However, when the aim is to make 
a precision measurement of the W mass (to better than one per mille), one cannot rashly 
neglect their possible contribution. Since we know of no "correct" procedure to calculate it 
the reasonable approach is to adopt a "good bet" default with a method to define a "band of 
uncertainty" . 



5.3.2 Flavour and kinematics specification 



It is natural to use the HEPEVT common block specification ||166|| to transfer flavour and kine- 
matics information from the electroweak generator to the QCD one. After all, the HEPEVT 
standard was devised specifically with this kind of tasks in mind. The original standard has 
been changed so that real variables are given in DOUBLE PRECISION. 

For the current interface, only NHEP, IDHEP and PHEP are actually mandatory. EW generator 
authors are invited to fill also the other information, such as mother-daughter pointers, but that 
is optional. Furthermore, any number of entries may be used in the event record to indicate the 
incoming e"'"e~ pair and intermediate states, but the only objects allowed to have status code 
ISTHEP= 1 are the two final fermion-antifermion pairs and an arbitrary number of photons. 
The fermions may be interspersed with photons in the listing, but the relative order of fermions 
is strict: 

1 one outgoing fermion, i.e. q/i'/ue] 

2 one outgoing antifermion, i.e. q/i^/Vf, 

3 another outgoing fermion; and 

4 another outgoing antifermion. 

The pairing of the outgoing fermions and antifermions should be done so that, when W'''W~ 
intermediate states can contribute, the pair 1 and 2 corresponds to a possible decay of the W"*", 
and the pair 3 and 4 to a possible decay of the W~. An example of an allowed order is uddu, 
while uudd is not correct. When a W'''W~ pair cannot contribute, the ordering should be 
instead made consistent with the decay of of one Z° to the pair 1 and 2, and another Z° to the 
pair 3 and 4, e.g. uucc. This way, coloured and uncoloured fermions can not be mixed in a 
pair, i.e. ue~I7ed is not an allowed ordering. 

Of course, adopting the fixed order above is not crucial, but it avoids the need for QCD 
generators to do a lot of rearrangements, and establishes a standard for the colour flow weights 
in the next section. 



5.3.3 Colour flow specification 



As was already mentioned above, the colour flow is not uniquely specified when both outgoing 
fermion pairs are of quark-antiquark type. A QCD event generator is therefore required to 
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make a choice. We propose the following procedure. 

The acceptance of a kinematical configuration by the electroweak generator (including spe- 
cific helicities for some generators) is based on the total squared amplitude, \A\'^, so this 
number is available "for free". Internally, a generator has access to the subamplitudes, A = 
Ai + A2 + . . . + Aji- Each subamplitude does correspond to a well-defined colour flow, so split 
the amplitudes into two classes, I and II, with I corresponding to colour singlets 1-1-2 and 3-1-4, 
and II to singlets 1-1-4 and 3-1-2. The total squared amplitude can then be written as 

\A\^ = \Ai + Auf = \Ai\^ + \Auf + 2Re(AiAi*i) = \Ai\^ + \An\^ + A . (21) 

This subdivision should be gauge invariant. 

Each electroweak event generator should return the three (positive) numbers \A\'^, lAi]"^ and 
I An p. Then the colour-suppressed interference term A is easily found as A = — |>lnp. 

A "good bet" approach to the colour assignment problem is for the QCD generator to 
neglect the interference term, and use the relative magnitude of \Ai\^ and |Anp to make a 
choice at random between the two possible colour flows. 

More sophisticated recipes are used for QCD processes like qg — > qg in HERWIG and 

Pythia, where the interference terms are split between the non-interference ones in accordance 
with the pole structure. However, such an approach presupposes a detailed study for each 
specific combination of allowed graphs, and so cannot be part of a generic interface. Should 
a generator provide such a subdivision, maybe as an option, it would be easy to represent by 
modified numbers \Ai\^ \Ai\^ + Aj so that A ^ A — Ai — An = 0. 

When A is nonvanishing, the uncertainty can be estimated by assigning the interference 
terms so that either class I or class II is maximized. Specifically, class I is maximized when the 
choice of colour flow is based on the relative magnitude of Ri and it!ii, where 

Ri =UiP + A Ru = lAiiP ifA>0 

(22) 

= \Ai\^ = max(0, \Au\^ + A) else 

and correspondingly with I-<->II for class II maximized. If the difference between these two 
extremes is small, then presumably the default procedure can be trusted. 

5.3.4 Further problems 

A number of potential problems exist, where the current approach may not be enough. These 
are discussed in the following. 

It has implicitly been assumed that the scale of perturbative QCD parton-shower evolution 
is set by the mass of the respective colour singlet. An example of a process where this need not 
be the case is e+e~ — > 7*/Z'^ — > uu — > uuZ° — > uudd. The uu pair here has a large original mass. 
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which is reduced by the emission of a Z . It is not clear whether the QCD radiation can be well 
approximated by that of the final uu mass, or whether the original mass is somehow felt e.g. 
by a larger rate of hard-gluon emission. A study of the matrix element for e''"e~ — > qq ^ qqZ°g 
would here be necessary. However, these graphs are not expected to give a major contribution, 
so presumably the uncertainty from this source is not significant. 

When QCD processes are introduced, interference terms need not be colour-suppressed. 
Specifically, the graph c+c" 7*/Z° — > uu uug uudd gives two colour singlets ud 
and du, just like a W^W^ intermediate state would. Therefore a suppression of interference 
contributions has to be based entirely on kinematical considerations. 

Since separation of quark and gluon jets is very difficult on an individual basis, also qqgg 
gives a background to four-fermion final states. Here, of course, there can be no interference 
with the other processes. 

The addition of parton showers to a QCD four-jet event, either qqgg or (qqg — ^)qqq'q', has to 
follow quite different rules from that of other four-fermion events, e.g. with respect to angular- 
ordering constraints in the parton shower. These rules have not yet been worked out for any of 
the QCD generators. The input that electroweak generators can give here is therefore not so 
meaningful. The main thrust in this area should be an improved matching between the matrix- 
element and the parton-shower strategies already present in QCD generators. Electroweak 
generators (if they contain QCD graphs) should therefore have the option of switching off all 
QCD contributions, i.e. (the amplitudes for) the graphs above. 

Some further input parameter may be required to specify whether QCD showers should be 
allowed also to involve the emission of photons. At LEP 1 we have learned that the "compe- 
tition" between photon and gluon emission is a not unimportant aspect, that tends to reduce 
the total amount of photon radiation compared to the no-QCD-radiation scenario. Something 
similar is likely to hold at LEP 2. However, the situation is far worse here, since the number 
of charged particles is much larger, and the presence of intermediate charged states (W'''W~) 
makes a subdivision of the full emission rate much more comphcated. One could therefore 
consider two extremes: 

• If an EW generator attempts to do the full job of photon radiation from all charged legs, 
then the QCD generator should not add further photon radiation. In fact, if anything, one 
may question whether the EW generator overestimated the amount of photon radiation 
off the quarks. 

• If an EW generator only claims to have initial-state photon radiation, then the QCD gen- 
erator could add final-state radiation inside each fermion-antifermion pair (also leptons, if 

implemented). This would still not be the full answer, but likely to be better than having 
no final-state radiation at all. (Since there is no unique, gauge-independent definition 
of final-state 7 radiation in four-fermion processes, the usefulness of such an approach 
should be checked from case to case.) 
Traditionally, QCD generators are not good at handhng the polarization of r's in the decay 
treatment. This is better done by dedicated r decay packages. Therefore EW generators that 
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do provide the spin of outgoing r's should give this information for the r entries in 

COMMON/HEPSPN/SHEP (4 , NMXHEP) 

using the standard conventions ||166|| . A flag could be set by the EW generator, and used by 
the QCD generator to inhibit it from decaying the r's. 

We remind the reader that the production vertices at the femtometer level may be of interest 
for physics such as colour reconnection and Bose-Einstein effects. If any generator should 
provide such output, the VHEP part of HEPEVT can be used to define vertices. The original 
objective was for vertices at the scale of mm, but also numbers of order 10~^^ mm could be 
stored with maintained precision so long as the primary event vertex is designed to be at the 
origin. 

5.3.5 Existing codes: 

Several interfaces now exist that are based on the philosophy outlined above. 

• Output from ExcALlBUR, with amplitude information for the different colour singlets. 
Can be obtained at 

pittp : //wwwcn. cern. ch/^charlton/excalibur/excalibur .html . 

• Input into HERWIG. Can be obtained at 

pittp : //suryall . cern. ch/users/seymour/herwig/ . 

• Input into Jetset. Can be obtained at 

pittp : / /thep . lu. se/tf 2/ staf f /torbj orn/test/main07 . f . 

• Input into Ariadne. Is part of the standard Ariadne distribution. 
Further information is available in the respective files. 

5.4 Systematic errors 

At LEP 2, several physics issues will involve hadronized quarks, both for QCD studies and for 
Electroweak measurements or searches. The following can be envisaged studies: 

• Establish the running of ttg from the Z pole to LEP 2 energies. 

• Measurement of hadronic cross-sections, e.g. e"^e^ ^ Z/7 ^ qq or e+e" — > W'''W~ 
Qi'l2l3Q4) discrimination between the two. 

• Reconstruction of jet-jet invariant masses in the above two processes, or even in e^e^ —>■ 
ZH qqbb. 

The first item seems the easiest case, as has been measured using a great variety of observables 
at the Z peak, with nearly infinite statistics. The variation with ^/s of well defined quantities 
such as energy-energy correlations or their asymmetry, or jet rates for a given ?/cut should be 
much less prone to systematic errors than their relationship to as itself. 
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Two difficulties can be expected here. First, the flavour composition of the sample will be 
different at LEP 2. In particular, the rate of bb production will decrease from 22% down to less 
than 10%. The fact that the specific fragmentation parameters for b quarks [ |17U|| have been 
measured at LEP 1 should be of great help. In order to extrapolate to higher energies, these 
results have to be incorporated in the simulation in one way or another. The three-jet rate has 



been used at the Z peak to test the universality of the strong coupling ||171|| with an accuracy of 
about 0.005. The argument can be turned around as, the sensitivity of a determination of as to 
flavour composition, leading to a rough uncertainty estimate of about 0.0005 on the difference 
in as from the Z peak to LEP 2. 

The second difficulty will arise when one tries to go from establishing the running to more 
quantitative estimates of it. The running will be compared to the expectation from the QCD 
fragmentation models. Given that the most popular generators are presently based on 0{as) 
exponentiated showers, one can rightfully challenge their capability to predict the -y/s evolution, 
because of missing higher orders. The solution to this issue will probably have to come from a 
better mapping of the shower models to second order matrix elements. 

The impact on acceptance corrections was limited at the Z peak by two positive factors: large 
statistics and limited initial state radiation. A simple event rotation technique ||172 was suf- 



ficient to reduce the uncertainty on event selection down to 10^^ or better. The precision 
required at LEP 2 for such studies is less stringent, statistical errors being at the level of 1%, 
so that the same method applied to high energy annihilation events should be adequate. One 
difficulty will arise from initial state radiation (ISR): the optimum sensitivity for electroweak 
effects is obtained by removing the radiative return to the Z peak using an s' cut. Most of the 
ISR photons being emitted at small angles, the invariant mass of the hadronic system has to 
be used to implement this cut. The issue here is to understand how accurately one can recon- 
struct an invariant mass from a system of boosted jets. An important experimental constraint 
can presumably be placed by using e+e~ — > Z + 7 —>■ qq-f7 events. However the issue of 
flavour dependence will come up again here, as the mass of the b quark and missing energy 
from neutrinos are expected to have sizeable effects on the jet angles and energies after a boost. 
A similar problem will be encountered when reconstructing W— qq invariant mass, where the 
difference in flavour composition is even more drastic. 

Finally one last but important issue is the discrimination between e"'"e~ Z/7 qq 
and e"'"e~ W"'"W~ ^1(12*13^4 events for the determination of the W+W~ cross-section at 
threshold [ |173[ . There is a finite probability that a four-jet event from the ffist process with 
two hard QCD-radiated partons will mimic the second process. In the present state of QCD 
generators with only 0{as) exponentiated showers, it is not obvious that the Monte Carlo gives 
the right answer. One way to obtain direct experimental information is to see how often a 
hadronic Z decay can be reconstructed as two heavy systems of 45 GeV mass and compare 
with the predictions of the fragmentation model. The extrapolation to the appropriate center- 
of-momentum energies and invariant mass requires a fragmentation Model. 
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In most of these problems, an experimental constraint can be found in e.g. Z decays. How- 
ever every time fragmentation event generators are needed to perform the necessary extrapo- 
lations. Evaluating the corresponding systematic errors has been performed traditionally by 
either i) varying some (well chosen) input parameters within "reasonable limits" or ii) compar- 
ing the results obtained when using two different models. Recently, a more complex situation 
has emerged for the analysis of the jet charge asymmetries in Z decays ||174|| . This is a clear 



example of an electroweak measurement performed using jets. The jet charge separations are 
ultimately obtained from a fragmentation model, upon which many constraints are imposed: 
measured production spectra for pions, kaons and baryons (p and A), resonances such as p, K* 
and rj, average jet charge measured from opposite hemisphere charge correlation, etc. Impos- 
ing these constraints immediately leads to extremely strong correlations among fragmentation 
parameters. In Jetset, it is possible to find enough parameters to describe very completely 
the production of each particle species. The weak points remain the transverse momentum 
distributions and the baryon spectra. In HERWIG fewer parameters are available and the is 
worse. Nevertheless the value of the electroweak asymmetry can be extracted for both models, 
with systematic errors related to the goodness of fit. A consistency check is supplied by the 
agreement of the values obtained from the two models within the systematics pertaining to 
each model. Similar procedures can be envisaged for measurements of electroweak quantities 
at LEP 2. 

To conclude, there is no unique prescription for evaluating systematic errors. In each prob- 
lem specific sources of errors and the corresponding fragmentation model parameters have to 
be found. Incorporating experimental constraints generally leads to very strong correlations 
among parameters, but this can be solved by for example using a combined linearized fit. The 
most general problem in extrapolating results obtained at the Z pole to LEP 2 will be the 
change in flavour composition. A mapping of the parton shower models to the second order 
matrix elements would be most useful. 
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6 Summary and Recommendations 



It is useful to remember the last words of the LEP 1 QCD generators report [Q: Due to the 
large uncertainty present in any realistic Monte Carlo, physics studies must he based on the use 
of at least two complete and independent programs. Nothing has been changed in this regard; 
QCD is still not solved and the need for models is as large as ever. 

The QCD generators of today may be considered more mature than the pre-LEP ones, 
in that they have successfully survived a number of experimental tests. However, there is 
always the danger that "incorrect" models do not just fade away — they are only modified and 
retuned for agreement. The increased energy lever arm provided by LEP 2 could give additional 
discrimination power, or at lease necessitate further fine tuning of programs. 

Furthermore, in comparing with the LEP 1 data, we see that no generator is perfect. De- 
pending on the physics area studied, it is therefore important to beware of generators with 
known shortcomings in that area. These shortcomings may indicate basic problems in the 
models, but could also come from further effects (e.g. higher-order matrix-element corrections) 
that authors never claimed to include. Generator authors are encouranged to sort out known 
problems in the light of LEP 1 experience, and in particular those with implications for LEP 2 
studies. For some areas, such as colour reconnection and Bose-Einstein effects, the modelling 
is only in its infancy, and further efforts obviously are required. 

The World Wide Web offers new opportunities to make programs accessible, including 
manuals, update notes, sample main programs and so on. To the extent authors did not yet 
adapt their distribution practices to the new opportunities, they are encouraged to do so. A 
common practice of having a "home page" for each generator will allow the construction of 
useful generator directories. 

Standardization is as important as ever, in order to avoid confusion among experimentalists 
required to run a multitude of different codes. We have here proposed modifications to the 
/HEPEVT/ standard, extensions and a few corrections to the PDG particle code, a standardized 
decay table and an interface between electroweak four-fermion generators and QCD generators. 
A continued dialogue about possible standards would be very useful. 
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